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reflection).  On  the  return  path,  the  aberrations  in  the  ampli¬ 
fier  are  compensated  for  by  the  conjugate  beam  produced  by  the 
SBS  conjugator.  Ideally,  the  result  is  an  amplified,  unaberrated 
output  beam. 

Even  though  the  SBS  conjugator  can  produce  a  high  quality 
conjugate  wavefront,  there  are  limitations  on  the  compensation  of 
the  RA:  (1)  possible  temporal  limitations  due  to  the  finite 

response  time  of  the  SBS  conjugation  process  (<10  ns), 

(2)  spatial  limitations  due  to  diffraction  and  clipping  in  the 
RA,  and  (3)  inability  to  compensate  for  aberrations  that  are 
different  on  the  initial  and  return  passes  through  the  RA  (nonre¬ 
ciprocal  aberrations) .  (Long  pulse  issues  characteristic  of 
excimer  systems  were  not  investigated  on  this  program.) 

Fluctuations  faster  than  100  MHz  appear  unlikely  in  RA 
systems  and  were  not  investigated  during  this  program.  The 
limits  on  correction  of  small-scale  fluctuations  from  diffraction 
and  clipping  have  been  extensively  evaluated  at  Hughes  Research 
Laboratories  for  both  amplifiers  and  oscillators.1’2’3  From 
these  studies  we  conclude  that  SBS  provides  excellent  correction 
for  aberrations  whose  rms  phase  is  <2 ir  rad,  provided  that  the 
Fresnel  number  (radius  squared  divided  by  length  times  wave¬ 
length)  is  >1 .  In  a  well-designed  system,  aberrations  from  dif¬ 
fraction  and  clipping  will  be  minimized;  their  effects  were  not 
reinvestigated  during  this  program. 

The  goal  of  this  program  was  to  determine  the  physical  pro¬ 
cesses  that  limit  the  compensation  of  RA  phase  aberration  by  SBS 
conjugators.  If  all  RA  aberrations  are  reciprocal  and  if  the 
conjugation  fidelity  of  the  SBS  process  is  sufficiently  good,  the 
wavefront  of  the  amplified  output  should  be  very  nearly  the  phase 
conjugate  of  the  input  reference.  Departures  from  this  ideal 
performance  are  caused  by  nonreciprocal  aberrations  and/or  imper¬ 
fect  performance  of  the  SBS  conjugator.  Nonreciprocal  aberra¬ 
tions  include  direct  nonreciprocal  phase  aberrations,  as  well  as 
all  amplitude  aberrations  due  to  nonuniform  gain-. 
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The  purpose  of  this  work  was  to  quantify  the  distortions 
caused  by  nonreciprocal  aberrations  in  the  RA  process  for  uniform 
and  distorted  pump  and  seed  waves.  Deleterious  effects  arise 
from  poor  replication  of  either  the  amplitude  or  phase  of  the 
input  seed.  Poor  phase  replication  leads  directly  to  phase 
errors  in  the  output,  while  poor  amplitude  replication  leads  to 
phase  errors  only  after  propagation.  It  is  the  coupling  of 
amplitude  and  phase  aberration,  by  diffraction,  that  leads  to  the 
varity  of  possible  effects.  Nonuniform  gain  caused  by  pump  aber¬ 
ration  and  pump  saturation  gives  rise  to  a  distorted  amplitude 
profile.  Because  the  gain  in  the  RA  is  distributed,  the  effect 
depends  on  the  scale  size  of  the  gain  nonuniformity  relative  to 
diffraction  (pump  and/or  seed) .  When  the  scale  size  of  the  gain 
nonuniformity  is  very  large  compared  with  diffraction  (uniform 
gain),  the  RA/SBS  rystem  works  very  well.  As  the  gain  nonuni¬ 
formity  scale  size  decreases,  its  effect  becomes  much  more  impor¬ 
tant.  In  the  limit  of  very  short  scale  size,  the  output  beam 
fidelity  improves  because  diffraction  tends  to  average  over  the 
gain  nonuniformity. 

A  perfect  conjugator  system  can  compensate  for  phase  aberra¬ 
tions  that  are  reciprocal,  such  as  thermal  blooming  and  refrac¬ 
tion  by  an  inhomogeneous  vibrational  excitation  density  (RIVED) . 
These  aberrations  are  correctible  if  the  resulting  changes  in 
refractive  index  occur  slowly  compared  to  the  round  trip  transit 
time  of  the  wavefronts.  The  conjugator  can  never  be  expected  to 
compensate  for  nonreciprocal  aberrations  (i.e.,  aberrations  such 
as  gain  saturation) . 

The  concern  addressed  under  this  program  was  that  nonrecip¬ 
rocal  aberrations  could  impair  the  ability  of  the  conjugator  to 
compensate  for  the  reciprocal  aberrations.  We  examined  this 
question  by  comparing  beam  fidelity  in  the  presence  of  nonrecip¬ 
rocal  aberrations  with,  and  without,  the  addition  of  fixed  recip¬ 
rocal  aberrations.  Fixed  reciprocal  aberrations  will  always 
exist  within  the  amplifier  as  a  result  of  thermal  effects  and  of 


distortions  in  optical  components.  These  aberrations,  however, 
are  difficult  to  quantify  and  are  generally  too  small  to  produce 
measurable  effects  in  our  experiments.  To  produce  large  effects 
and  to  demonstrate  the  ability  of  the  conjugator  to  compensate 
for  large  reciprocal  aberrations,  we  included  a  fixed  aberration 
plate  in  front  of  the  amplifier.  This  plate  is  used  to  simulate 
atmospheric  effects  as  well  as  aberrations  in  external  optics  and 
reciprocal  aberrationj  in  the  amplifier. 

To  quantify  the  results  of  this  study,  we  compare  the  beam 
quality  of  the  aberrated  system  to  an  identical  system  where  the 
aberrator  is  removed.  This  second  unaberrated  system  measures 
the  effect  of  the  nonreciprocal  phase  errors;  comparison  with  the 
aberrated  system  allows  evaluation  of  the  ability  of  the  conjuga¬ 
tor  to  compensate  for  the  reciprocal  errors. 

The  results  of  this  program  can  be  summarized  as  follows: 

•  Identified  five  sources  of  nonreciprocal  aberrations  in 
RA,  namely,  dispersion,  transverse  mode  four-wave  mixing 
(pump  replication) ,  longitudinal  mode  four-wave  mixing, 
pump  saturation ,  and  RIVED . 

•  Estimated  the  conditions  under  which  each  effect  can 
become  important. 

•  Modified  an  existing  propagation  code  to  include  the 
effects  of  nonreciprocal  aberrations  and  used  the  code 
to  model  the  RA/SBS  conjugator  system. 

•  Using  the  code,  we  found  that  the  conjugator  can  compen¬ 
sate  for  reciprocal  phase  errors  in  almost  all 
situations . 

•  A  few  situations  were  found  where  the  aberrated  system 
beam  fidelity  droped  by  a  factor  of  ~3  from  that  of  the 
unaberrated  system.  In  these  cases,  however,  we  were 
modeling  such  highly  aberrated  systems  that  the  beam 
fidelities  would  have  been  ~3  orders  of  magnitude  lower 
if  the  system  had  no  conjugator. 

•  Experimental  investigations  were  carried  out  using  a 
hydrogen  Raman  amplifier  and  the  second  harmonic  of  a 
Nd:YAG  laser  as  the  pump  beam.  These  experiments  were 
carried  out  with  both  single-mode  and  highly  multimode 
laser  sources . 
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•  Using  both  the  multimode  and  single-mode  laser  sour 
we  observed  a  steady  decrease  in  conjugate  beam  fid 
with  increasing  pump  energy.  We  believe  this  resul 
from  nonreciprocal  aberrations  due  to  nonuniform  ga 
and  pump  saturation. 

•  Single-mode  laser  results  were  compared  to  estimate 
from  the  numerical  simulation.  The  two  results  wer 
found  to  be  consistent. 

This  report  consists  of  five  sections.  The  first  sectio 
provides  a  brief  introduction  and  description  of  the  program 
Section  2  contains  an  analytical  task  aimed  at  determining  t 
processes  leading  to  nonreciprocal  aberrations  in  RAs ,  their 
relative  strengths,  and  the  means  to  minimize  their  effects. 
Section  3  describes  numerical  modeling  performed  using  a  thr 
dimensional  (3-D)  wave  optics  simulation  on  the  double-pass 
RA/SBS  system.  These  results  show  the  behavior  of  these  sys 
under  a  variety  of  operating  conditions,  in  addition  to  prov 
results  for  comparison  with  experiments.  Experiments  to  dem 
strate  the  nonreciprocal  phase  errors  and  to  compare  beam  fi 
ties  with  the  numerical  model  are  described  in  Section  4.  I 
Section  5,  we  compare  the  numerical  and  experimental  results 
draw  conclusions,  and  indicate  directions  for  future  work. 
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SECTION  2 


THEORY  OF  PHASE  CONJUGATION  AND  APPLICATION  TO  RAMAN  AMPLIFIERS 

2.1  THEORY  OF  RAMAN  AMPLIFIERS  AND  SBS  PHASE  CONJUGATION 

Ordinary  Raman  scattering  is  a  useful  tool  for  studying  the 
vibrational  energy  levels  of  molecules  and  optical-branch  lattice 
vibrations  in  crystals.  At  high  optical  power,  the  scattering 
process  becomes  stimulated  and  a  significant  portion  of  the  inci¬ 
dent  radiation  can  be  scattered.  Spectral  analysis  of  the  scat¬ 
tered  radiation  reveals  the  existence  of  frequencies  that  are 
shifted  down  from  the  incident  frequency  by  amounts  equal  to  the 
vibrational  frequencies  of  the  material  being  irradiated.  This 
scattering  process  is  referred  to  as  "Stokes  scattering."  The 
scattering  event  is  illustrated  in  Figure  2.  In  the  figure,  the 
molecule  is  initially  in  the  ground  state  (is  =  0) .  An  incident 
photon  at  frequency  wp  is  absorbed,  and  a  Stokes  photon  is  emit¬ 
ted  at  a  lower  frequency  w  s  .  The  frequency  mismatch  uu  =  wp  -  uis 
is  equal  to  the  vibrational  frequency  of  the  molecule.  To  con¬ 
serve  energy,  the  molecule  is  transferred  to  an  excited 
vibrational  level  (is  =  1)  . 
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Figure  2.  The  Stokes  scattering  process,  in  which  an 

input  pump  photon  at  frequency  u»p  is  absorbed 
and  a  Stokes  photon  at  frequency  ws  =  ui  p  -  is 
emitted.  The  excess  energy  is  taken  up  by 
exciting  the  molecule  to  its  first  vibrational 
level  (is  -  1)  . 
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Our  concern  in  this  section  is  the  process  by  which  energy  is 
transferred  from  pump  photons  to  Stokes  photons,  and  the  way  in 
which  this  process  affects  the  fidelity  of  the  Stokes  (seed) 
wavefront.  The  effects  of  interest  to  this  program  can  be  des¬ 
cribed  using  the  heuristic  model  of  Raman  scattering.  In  this 
model,  the  induced  dipole  moment  of  the  molecule  is  given  by: 

/*  =  e0ctE  ,  (1) 

where  E  is  the  applied  electric  field  and  a  is  the  polarizability 
of  the  molecule.  Clearly,  the  polarizability  of  the  molecule 
depends  on  the  intramolecular  spacing  and  therefore  will  be  modu¬ 
lated  as  the  molecule  vibrates.  Assuming  that  a  single  vibra¬ 
tional  mode  of  the  molecule  is  important,  the  dipole  moment  can 
be  expressed  as 


where  Q  is  the  vibrational  coordinate  of  the  dominant  mode.  The 
equation  of  motion  for  the  normal  mode  coordinate  Q  is 


where  Mj/  is  the  effective  mass  of  the  vibrational  mode,  7^  is  a 
phenomological  damping  constant,  and  u/y  is  the  normal  oscillation 
frequency  of  the  mode.  E2  in  the  equation  above  contains  terms 
that  oscillate  at  the  sum  and  difference  frequencies  of  the  pump 
and  Stokes  waves.  Only  the  difference  frequency  terms  oscillate 
near  the  normal  mode  frequency  and  are  able  to  drive  the  normal 
mode  vibration.  The  amplitude  of  this  vibration  is  given  by: 


Qi/  = 


i£°  raa) 

IdQj 


(4) 


where  Qj,,  Ep,  and  Es  are  the  Fourier  coefficients  of  the  normal 
modes  and  the  applied  fields  [i.e.,  Q  =  Q^exp (-iw^t)  + 

Q^exp  (iu/^t)  ,  etc.].  As  discussed  above,  this  vibrational  mode 
drives  the  molecular  polarizability  which,  as  a  result  of  the 
applied  electric  field,  drives  the  induced  dipole  moment  of  the 
molecule.  Using  the  result  of  Eqs .  (2)  and  (4),  we  find  dipole 

moments  oscillating  at  both  the  Stokes  and  pump  frequencies. 

This  gives  rise  to  a  macroscopic  polarization  density  at  the 
Stokes  frequency  given  by: 

2 

Pe  o  (da')  2 

p.  -  III)  m;  e.  • 

This  macroscopic  polarization  P,  drives  the  Stokes  wave,  giving 
rise  to  the  optical  gain  in  the  Raman  amplification  process. 

From  Eq.  (5)  the  following  conclusions  can  be  drawn: 

•  The  Raman  gain  depends  only  on  the  pump  intensity 
(la  I E „ I  2 )  ,  and  not  its  phase.  It  follows  that  pump 
phase  aberration  has  no  effect  unless  the  effects  of 
diffraction  cause  the  phase  aberration  to  manifest  an 
intensity  aberration. 

•  The  Raman-driven  polarization  (Pt)  has  the  spatial  pro¬ 
file  of  the  seed  wave  (i.e. ,  P*  a  Es) .  Therefore,  the 
amplified  seed  wave  can  be  a  high  fidelity  replication 
of  the  input  seed  wave  if  the  pump  profile  is  flat 

(epe;>. 

In  this  manner,  as  is  known,  the  Raman  process  can  be  made  to 
transfer  energy  from  an  aberrated  (phase)  pump  beam  to  an  unaber¬ 
rated  seed  beam. 

This  simple  model  predicts  high  fidelity  amplification.  In 
real  systems,  however,  other  effects  must  be  taken  into  account. 
These  include  diffraction,  pump  depletion  (gain  saturation), 
saturation  of  the  Raman  medium,  and  dispersion.  Diffraction 
leads  to  intensity  aberrations  on  the  pump  beam  and  therefore  to 
nonuniform  gain,  giving  rise  to  a  poorly  replicated  seed  beam. 


In  regions  where  the  seed  beam  is  intense,  the  pump  becomes  more 
highly  depleted,  also  giving  rise  to  a  nonuniform  gain.  This  is 
especially  important  for  high  gain  and  conversion  efficiency.  In 
addition  to  gain  saturation,  the  medium  can  be  saturated,  which 
occurs  if  a  significant  fraction  of  the  molecules  are  left  in 
their  excited  state  (v  =  1) .  Finally,  dispersion  effects  from 
nonuniform  gain  and/or  saturation  of  the  medium  can  distort  the 
amplified  seed  beam.  In  reality,  all  of  these  effects  are  pre¬ 
sent  simultaneously,  and  it  is  their  mutual  interaction  that 
distorts  the  amplified  seed.  The  purpose  of  this  program  is  to 
determine  the  important  effects,  the  conditions  under  which  they 
seriously  degrade  the  amplified  beam,  and  the  extent  to  which 
phase  conjugation  can  improve  the  output  beam  quality. 

SBS  production  of  a  phase  conjugate  wave  can  be  understood  by 
considering  three  properties  of  SBS: 

•  SBS  grows  from  noise. 

•  The  total  SBS  gain  is  large  (>e2s). 

•  The  local  gain  coefficient  in  any  thin  layer  of  the 
medium  is  proportional  to  the  pump  intensity. 

Starting  at  an  infinitesimal  layer  at  the  back  of  the  SBS  cell, 
the  noise  radiation  that  matches  the  pump  intensity  undergoes  the 
largest  gain.  Of  course,  this  radiation  is  not  phase  conjugate 
to  the  pump  and  does  not  have  a  coherent  phase  front.  However, 
if  the  SBS  interaction  length  is  long  enough  and  if  the  decor¬ 
relation  length  of  the  pump  intensity  in  the  direction  of  propa¬ 
gation  is  small  enough,  this  radiation  diffracts  as  it  propagates 
to  a  new  infinitesimal  layer  closer  to  the  entrance  of  the  cell. 
At  this  layer,  the  amplified  noise  undergoes  another  amplifica¬ 
tion/filtration  process.  Only  the  portion  of  the  noise  with  the 
proper  phase  front  to  diffract  into  the  pump  intensity  profile  in 
the  next  layer  will  be  amplified.  As  this  process  continues 
through  a  large  interaction  length,  the  phase  conjugate  field  has 
a  gain  coefficient  approximately  twice  that  of  the  nonconjugate 
component,  dominating  the  nonconjugate  component. 


2.2  SOURCES  OF  NONRECIPROCAL  ABERRATIONS 


At  least  five  sources  of  nonreciprocal  aberrations  exist  in 
RA  systems,  including  dispersion,  four-wave  mixing,  pump  replica¬ 
tion,  amplifier  saturation,  and  self-focusing.  These  mechanisms 
depend  on  details  of  the  pump,  the  Stokes  seed,  and  the  amplifi¬ 
cation  geometry.  While  these  mechanisms  have  been  studied  in 
past  work,  the  ability  of  a  double-pass  conjugation  system  to 
compensate  for  their  aberrations  has  not. 

Dispersion  in  a  RA  is  analogous  to  anomalous  dispersion  near 
a  two-level  gain  or  absorption  line,  and  can  lead  to  phase  per¬ 
turbations  if  the  pump  and  Stokes  are  not  tuned  exactly  to  the 
Raman  transition.4,5’6  If  the  pump  intensity  is  uniform,  the 
effect  gives  rise  to  a  uniform  shift  of  the  seed  phase  and  does 
not  reduce  the  beam  fidelity.  A  nonuniform  pump  beam  introduces 
a  nonuniform  phase  shift,  giving  rise  to  a  direct  phase  aberra¬ 
tion  of  the  seed  beam.  This  aberration  occurs  for  multimode  pump 
and  Stokes  seed,  or  a  single-mode  Stokes  seed  that  is  off-line 
center  (because  of  incomplete  compensation  of  the  SBS  shift) . 

The  magnitude  of  the  effect  can  be  evaluated  from  the  Kramers- 
Kronig  relationship,  which  shows  that  the  gain  and  the  refractive 
index  are  intimately  related.  As  shown  in  Figure  3,  the  refrac¬ 
tive  index  is  decreased  on  the  low  frequency  side  of  the  gain  and 
increased  on  the  high  frequency  side.  The  phase  aberration  pro¬ 
duced  depends  on  the  spatial  nonuniformity  of  the  Raman  gain 
coefficient,  which  is  proportional  to  the  pump  intensity. 

Because  the  pump  and  seed  waves  are  crosspolar i zed  on  the  first 
pass  through  the  amplifier,  little  Raman  gain  is  experienced,  as 
based  on  a  calculation7  of  H2  polarizability  as  a  function  of 
internuclear  separation.  From  this,  we  estimate  that  the  cross- 
polarized  Raman  gain  is  <0.1  of  the  copolarized  gain.  Hence,  the 
nonlinear  dispersion  aberration  is  nonreciprocal,  occurring  only 
during  the  second  pass  through  the  amplifier.  If,  in  addition, 
the  Raman  transition  saturates,  the  phase  shift  from  dispersion 
also  saturates,  depending  on  the  pump  intensity  'in  a  highly  non¬ 
linear  manner.8 


Figure  3.  Plot  of  the  index  perturbation  around  a 
typical  gain  line. 

The  maximum  index  perturbation  occurs  where  the  gain  has 
fallen  off  by  a  factor  of  2.  At  this  point,  the  index  perturba¬ 
tion  is  ~(gIL)X/2L,  where  gIL  is  the  peak  small  signal  Raman  gain 
coefficient.  It  follows  that  the  maximum  integrated  phase  error 
for  propagation  across  the  amplifier  is  ~7r  (gIL)  .  At  typical 
gains,  several  waves  of  phase  shift  could  result;  however,  it  is 
only  the  spatial  nonuniformity  of  the  phase  shift  (gain)  that 
leads  to  aberration.  In  cases  where  this  dispersion  effect  would 
be  significant,  the  effect  of  the  nonuniform  gain  on  the  inten¬ 
sity  profile  is  much  more  severe,  masking  that  from  dispersion- 
Because  of  its  nonreciprocal  nature,  this  dispersion  effect  can¬ 
not  be  compensated  by  the  phase  conjugate  mirror.  The  effect  of 
this  in  the  RA/SBS  conjugator  system  is,  however,  not  clear. 

Four-wave  mixing  of  multiple  spatial  (transverse)  modes  or 
multiple  longitudinal  modes10  can  lead  to  nonreciprocal  aberra¬ 
tions  on  the  Stokes  beam.  Consider,  first,  four-wave  mixing  of 
spatial  modes.  This  mixing  process  is  generally  described  using 
the  nonlinear  susceptibility  tensor  .  The  nonlinear  polari¬ 

zation  oscillating,  at  the  Stokes  frequency  is  given  by: 


Mm  ZZZ&Ztm 


P*  2  =  6X<3>EplE;2Eslexp[(kpl-kp2+k#1-k,2)*x] 


(6) 


where  the  subscripts  s  and  p  represent  Stokes  and  pump,  respec¬ 
tively.  When  a  large  number  of  spatial  modes  (off-axis  k  vec¬ 
tors)  are  present,  phase  matching  of  the  pump  and  Stokes  can 
occur  in  two  ways:  (1)  A  pump  mode  (Epl)  combines  with  its  com¬ 
plex  conjugate  (Ep2  =  Epl,  kp2  =  kpl)  and  a  Stokes  mode  to 
produce  gain  for  that  mode  (k#2  =  ksl),  as  shown  by  the  wave- 
vector  diagram  in  Figure  4(a) .  When  pump  modes  combine  in  this 
manner,  no  aberration  is  produced.  (2)  Two  different  pump  modes 
(Ep2  ^  Epl,  kp2  ^  kpl)  combine  to  generate  an  off-axis  Stokes 
mode  (ks2  ^  ksl),  as  shown  in  Figure  4(b).  With  this  latter 
occurrence,  the  one  Stokes  mode  is  mixed  with  another,  giving 
rise  to  an  aberration  in  (..he  intensity  and  phase  of  the  seed 
wave.  From  a  simple  analysis11  that  neglects  depletion  of  the 
pump  wave,  transverse  mode  four-wave  mixing  is  expected  to 
significantly  aberrate  the  Stokes  wave  if 

2*(gIL)F  «  (XDL) 2  <<  e»lL  ,  (7 
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Figure  4.  Four-wave  mixing  with  multi-transverse 
mode  pump  and  Stokes  waves . 
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where  gIL  is  the  small  signal  (log)  gain  in  the  RA,  F  is  the 
Fresnel  number  of  the  pump  beam  (F  =  D2/4XL,  with  D  the  pump  beam 
diameter,  X  the  wavelength,  L  the  interaction  length),  and  XDL  is 
the  times  diffraction  limited  (XDL  -  D/a,  where  a  is  the  trans¬ 
verse  coherence  length,  or  XDL  -  DA0/X,  where  L6  is  the  beam 
divergence) .  The  lower  limit  on  XDL  simply  states  that  an 
unaberrated  pump  beam  will  not  aberrate  the  Stokes  wave.  The 
upper  limit  expresses  the  fact  that  the  inhomogenities  of  a  suf¬ 
ficiently  aberrated  pump  beam  will  "average  out." 

If  the  pump  and  Stokes  consist  of  a  single  spatial  mode  but 
multiple  longitudinal  modes,  then  the  nth  pump  and  pth  Stokes 
modes  can  combine  with  the  mth  pump  mode  to  amplify  the  qfc h 
Stokes  mode.  If  the  nth  and  mth  pump  modes  have  different 
phases,  the  phase  difference  A0nm  =  <p n  -  <pm  is  transferred  to  the 


ql h  Stokes  mode  and  an  aberration  is  produced.  If  the  pump  phase 


difference  (A0nm)  is  spatially  uniform,  there  is  no  reduction  in 
beam  fidelity.  If,  however,  the  phase  difference  depends  on  the 
transverse  direction,  an  aberration  is  produced.12’13  It  has 
been  suggested  that  <p  n  and  $m  are  likely  to  be  the  same  for 
multi-longitudinal  mode  radiation  from  excimer  lasers14  and  that 
therefore  four-wave  mixing  of  longitudinal  modes  may  not  be  an 
important  source  of  aberrations . 

Pump  replication  is  another  name  for  transverse  mode  four- 
wave  mixing,  referring  to  a  process  by  which  the  amplitude  and 
phase  of  the  Stokes  field  tend  to  replicate  the  pump  amplitude 
and  phase.  Replication  results  because  the  Stokes  gain  is  pro¬ 
portional  to  the  pump  amplitude,  and  hence  the  Stokes  wave 
becomes  intense  in  the  same  regions  that  the  pump  is  intense.  If 
this  process  is  accompanied  by  diffraction,  the  pump  amplitude 
and  phase  can  be  transferred  to  the  Stokes  wave . 1 5 ’ 1 ® ’ 1 7  Pump 
replication  for  the  forward  stimulated  scattering  processes  is 
analogous  to  phase  conjugation  in  backward  stimulated  scattering. 
This  means  the  gain  coefficient  for  replication  of  the  pump 
intensity  and  phase  is  about  a  factor  of  2  larger  than  the  gain 
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coefficient  of  a  plane  or  image-bearing  wave.  This  process  is 
one  of  the  ultimate  limits  on  the  injection  level  of  a  RA.  For 
example,  if  the  amplifier  gain  is  103  for  the  seed  wave,  the  gain 
for  noise  (which  matches  the  pump  amplitude  and  phase)  could 
approach  10®,  and  noise  at  levels  of  10' 3  of  the  seed  competes 
with  the  seed. 

While  a  double-pass  phase  conjugation  system  cannot  compen¬ 
sate  for  pump  replication,  it  can,  in  an  alternative  system, 
eliminate  the  problem.  This  objective  can  be  accomplished  by 
finding  a  means  to  introduce  aberrators  in  the  amplifier  that 
affect  only  the  Stokes  wavelength,  or  that  are  highly  dispersive 
so  they  affect  the  Stokes  wave  and  the  pump  wave  quite 
differently.  The  effect  of  these  aberrators  on  the  injected  seed 
is  completely  compensated  by  the  conjugator  (assuming  that  the 
beams  are  not  clipped  by  the  system  aperature) .  The  pump  and 
replication  waves,  however,  remain  uncompensated.  It  follows 
that  the  gain  length  for  the  replication  mode  is  limited  by  the 
aberrator  spacing  rather  than  the  amplifier  length.  Because  the 
gain  length  appears  in  the  exponent,  small  reductions  in  the  gain 
length  can  lead  to  large  reductions  in  the  replication  wave  gain. 

A  related  effect  is  pump  saturation,18  which  limits  the  gain 
in  regions  where  most  of  the  pump  energy  has  already  been  con¬ 
verted  to  the  Stokes  wavelength.  As  a  result,  the  amplitude 
profile  of  the  Stokes  wave  tends  to  match  the  pump,  rather  than 
the  amplitude  of  the  seed  wave.  Subsequent  propagation  converts 
this  amplitude  aberration  into  a  phase  aberration.  Pump  satura¬ 
tion  differs  from  pump  replication  because  the  pump  energy  is 
depleted,  causing  the  Stokes  and  pump  profiles  to  match  at  a 
single  point  within  the  amplifier.  That  is,  the  saturation 
effect  is  localized  and  cannot  couple  with  diffraction  to  produce 
pump  replication. 

Refraction  by  an  inhomogeneous  vibrational  excitation  density 
(RIVED)  refers  to  a  process  by  which  the  Raman  amplification 
process  changes  the  Raman  active  medium.  We  consider  here  the 
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three  ways  the  system  is  affected  by  promoting  significant 
numbers  of  molecules  to  the  excited  vibrational  state: 

•  Promotion  of  molecules  to  the  excited  state  means  that 
they  are  no  longer  available  for  the  Raman  process, 
leading  to  a  saturation  of  the  gain  medium.  If  the 
saturation  is  nonuniform,  the  intensity  of  the  amplified 
seed  will  be  aberrated. 

•  Molecules  in  the  excited  vibrational  state  have  a  dif¬ 
ferent  polarizability  from  those  in  the  ground  state, 
leading  to  an  index  perturbation  and  therefore  a  phase 
aberration  of  the  seed. 

•  The  Raman  process  deposits  significant  vibrational 
energy  into  the  medium.  As  this  energy  thermalizes, 
strong  blooming  and  beam  breakup  can  occur. 

The  latter  two  lead  to  phase  aberrations  that  are  generally 
reciprocal.  Only  in  the  short  pulse  regime,  where  transit  time 
is  important,  will  the  effects  be  nonreciprocal. 

Saturation  of  the  gain  medium  is  nonreciprocal  and  can  lead 
to  large  intensity  aberrations.  In  a  hydrogen  Raman  cell  with  a 
gas  density  of  100  amagats ,  the  molecular  density  is  ~1021  cm*3. 
Assuming  a  typical  interaction  length  of  100  cm,  it  would  take 
~1023  photons/cm2  to  saturate  the  medium.  This  corresponds  to 
roughly  106  J/cm2  for  visible  radiation.  Clearly,  other  damage 
mechanisms  become  important  long  before  these  levels  are 
attained . 

Other  nonlinearities  in  the  Raman  cell  can  be  important 
sources  of  nonreciprocal  phase  errors.  These  include  self-focus¬ 
ing,19,20,21  higher  order  Stokes  generation,  and  other  self- 
induced  aberrations  of  the  seed  beam.  Self-focusing  is  the  pro¬ 
cess  by  which  an  optical  beam  induces  a  change  of  refractive 
index  in  the  surrounding  media,  altering  its  propagation  charac¬ 
teristics.  If  the  induced  index  change  is  positive,  the  beam 
"self-focuses."  Because  the  seed  beam  is  much  more  intense  after 
amplification,  the  process  is  nonreciprocal.  In  gases,  the 
dominate  nonlinear  mechanism  for  this  effect  is  electrostriction . 


For  hydrogen  at  ~100  amagats ,  the  critical  power  for  self-focus- 
ing  (sel f -trapping)  is  ~1  MW/cm2.  Fortunately,  the  response  time 
for  the  effect  is  limited  to  the  transverse  beam  dimension 
divided  by  the  phonon  velocity,  or  about  5  fis  for  our  experi¬ 
ments,  and  even  longer  in  high  Fresnel  number  excimer  laser 
systems.  Our  experimental  pulse  lengths  are  less  than  20  ns,  and 
sel f-f ocusing  should  not  be  important.  Higher  order  Stokes 
generation  is  only  important  for  high  intensity  Stokes  beams .  In 
fact,  it  is  only  important  for  Stokes  energies  approaching  or 
exceeding  the  pump  energy  (due  to  shorter  interaction  length  in 
well-designed  systems).  Higher  order  Stokes  generation  can  lead 
to  a  large  nonlinear  loss  for  the  seed  wave,  depleting  the  seed 
where  it  is  most  intense.  We  do  not,  however,  believe  it  could 
have  been  important  in  our  experiments. 

2.3  NONRECIPROCAL  PHASE  ERROR  DISCRIMINATORS 

It  is  the  mutual  interaction  of  diffraction  and  the  effects 
mentioned  above  that  gives  rise  to  nonreciprocal  aberrations  and 
distortion  of  the  seed  beam.  For  example,  pump  replication 
requires  that  phase  aberration  and  diffraction  of  the  pump, 
spatially  nonuniform  gain,  and  intensity  aberration  and  diffrac¬ 
tion  of  the  seed  all  occur  simultaneously.  Because  of  this 
mutual  interaction,  the  identification  of  a  single  effect  or 
process  is  difficult,  if  not  impossible.  In  numerical  simula¬ 
tions,  the  effects  described  above  can  be  inserted  independently; 
in  this  manner  such  processes  as  pump  replication  can  be  analyzed 
without  competition  from  such  other  processes  as  pump  depletion. 
In  experimental  investigations,  the  situation  is  much  more  com¬ 
plicated.  The  initial  conditions  (such  as  beam  profile)  cannot 
be  arbitrarily  set  and  the  various  material  parameters  cannot  be 
arbitrarily  adjusted.  To  identify  the  important  effects  and 
processes,  we  need  to  find  the  system  dependence  on  the  adjust¬ 
able  parameters.  These  include: 


•  Raman  gain . 

•  Dispersion. 

•  RIVED . 

•  Pump  intensity. 

•  Seed  intensity. 

•  Injection  ratio. 

•  Internal  phase  aberrators. 

•  External  phase  aberrators . 

•  Raman  medium. 

•  Polarization. 

•  Temporal  overlap. 

•  Compensation  of  SBS  frequency  shift. 

The  first  three  of  these  can  be  arbitrarily  set  only  during 
numerical  simulations.  The  effect  of  these  parameters  on  the 
system  performance  will  be  discussed  in  the  following  sections 
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SECTION  3 


NUMERICAL  SIMULATIONS 

For  this  program,  we  modified  an  existing  computer  simulation 
of  a  RA/adaptive  optics  configuration  in  which  the  adaptive 
optics  "look"  through  the  RA .  The  techniques  used  in  this  model¬ 
ing  are  described  in  detail  in  Appendix  A.  In  this  section  of 
the  report,  we  summarize  the  salient  features  of  the  modified 
program,  indicate  how  we  used  it  for  these  studies,  and  describe 
the  results . 

3.1  THE  COMPUTER  MODEL 

Our  computer  model  has  been  implemented  for  3-D  propagation; 
that  is,  there  is  the  propagation  direction  and  two  dimensions 
transverse  to  the  propagation  direction.  The  code  does  not 
inherently  include  time-dependent  phenomena.  We  have,  however, 
simulated  time-dependent  phenomena,  such  as  RIVED,  in  an 
approximate  manner.  The  simulation  solves  the  coupled  nonlinear 
differential  equations  satisfied  by  the  pump  and  first  Stokes 
wave.  A  converging  coordinate  formulation  that  allows  us  to 
model  both  focused  and  collimated  wave  propagation  is  used  in  the 
simulation . 

The  optics  configuration  of  our  model  is  shown  in  Figure  5. 

It  consists  of  a  RA,  a  phase  conjugate  mirror,  phase  aberration 
screens,  and  the  optical  beams .  The  incoming  Stokes  seed  is 
assumed  to  be  uniform  (flat  top  or  truncated  Gaussian,  as  dis¬ 
cussed  below) .  When  this  beam  passes  through  the  aberration 
screen,  it  acquires  an  aberrated  phase,  which,  depending  on  the 
degree  of  aberration,  can  manifest  an  intensity  aberration  caused 
by  diffraction.  The  phase  screen  is  used  to  model  possible 
atmospheric  aberrations  and,  more  importantly,  to  impose  an 
adjustable  reciprocal  aberration  on  the  seed  wave.  The  seed  wave 
is  linearly  propagated  through  the  RA  cell  to  the  phase  conjugate 
mirror.  After  conjugation  or  reflection,  the  Stokes  seed 
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Figure  5 


The  optical  configuration  used  in  the 
simulation  code  to  model  the  characteristics 
of  the  double-pass  RA/SBS  conjugator  system. 


retraces  its  path  through  the  RA,  this  time  undergoing  gain  (non¬ 
linear  propagation)  from  the  copropagating  pump  wave.  After 
retraversing  the  phase  aberration  screen,  the  Stokes  seed  is 
examined  in  the  far-field.  The  pump  beam,  incident  from  the 
bottom  of  the  figure,  passes  through  an  independent  aberration 
screen  before  entering  the  RA .  Because  of  diffraction,  the  pro¬ 
pagation  distance  between  the  pump  aberrator  and  the  RA  allows 
intensity  aberrations  to  develop  on  the  pump.  This  development 
is  important  because  the  Raman  gain  process  is  automatically 
phase  matched,  and  phase-only  aberrations  of  the  pump  do  not 
aberrate  the  seed  wave . 

The  phase  distortion  plates  are  obtained  by  a  discrete 
Fourier  transform  process.  We  use  the  fact  that  the  Fourier 
transform  of  the  square  root  of  the  power  spectrum,  multiplied  by 
a  delta-correlated  random  function,  yields  a  distribution  with 
the  desired  correlation  properties.  A  2-D  array  containing 
Gaussian  random  numbers  (with  zero  mean  and  unit  variance)  is 
multiplied  by  a  spatial  frequency  envelope  described  by  a  modi¬ 
fied  Von  Karman  spectrum  of  the  form 
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where  K  is  the  spatial  frequency,  K0  =  27r/L0,  L0  is  the  outer 
scale  of  the  turbulence,  and  LM  is  the  width  of  the  mesh.  This 
gives  a  wave  with  the  characteristics  of  atmospheric  turbulence 
in  the  Kolmogorov  model.22  We  use  an  outer  scale  of  twice  the 
beam  diameter.  Figure  6  shows  a  typical  phase  aberration  profile 
that  is  multiplied  by  the  desired  RMS  phase  error.  For  the  pump 
wave,  the  phase  screen  is  used  with  propagation  to  produce  an 
intensity  aberration  before  the  RA .  For  the  seed  wave,  the  phase 
screen  is  used  to  simulate  all  reciprocal  phase  errors  the  seed 


beam  encounters. 


These  errors  include  those  originating  from  atmospheric  tur¬ 
bulence  and  optics,  as  well  as  those  arising  from  reciprocal 
aberrations  within  the  Raman  cell.  More  importantly,  the  seed 
phase  screen  is  used  to  impose  an  adjustable  reciprocal  aberra¬ 
tion  to  demonstrate  the  ability  of  the  conjugator  to  compensate 
for  reciprocal  aberrations  in  the  RA  system. 


Figure  6. 


Pseudo-color  plot  of  a  typical  phase  screen 
used  to  distort  the  pump  and  Stokes  waves. 
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Linear  propagation  is  simulated  by  a  solution  of  the  propaga¬ 
tion  integral  that  uses  either  a  Hermi te-Gaussian  transform  or 
the  discrete  Fourier  transform  technique  described  in  Appendix  A. 
We  found  the  Fourier  transform  technique  to  be  more  efficient 
when  the  beam  is  highly  aberrated.  Propagation  is  accompanied  by 
a  transformation  that  adjusts  the  size  of  the  beam  mesh  so  that 
the  ratio  of  the  beam  size  to  the  mesh  size  stays  roughly  con¬ 
stant  as  the  beam  diffracts.  This  transformation  is  normalized 


to  the  diffraction  angle  of  a  Gaussian  beam.  A  convergence  para¬ 
meter  a  is  chosen  to  increase  the  rate  at  which  the  mesh  size 
changes  to  account  for  the  larger  diffraction  angle  of  our  dis¬ 
torted  beams.  Typically,  a  value  of  2  (2X  the  Gaussian  beam 
diffraction  angle)  is  needed  to  allow  a  badly  aberrated  seed  to 
reach  the  conjugator  without  clipping  the  mesh  edge. 

Analytically,  a  perfect  conjugator  is  easy  to  implement;  it 
is,  however,  not  clear  that  the  SBS  process  always  yields  perfect 
conjugation.  For  the  purposes  of  this  analysis,  we  assumed  the 
conjugator  operated  in  one  of  three  modes: 


•  Perfect  phase  and  amplitude  conjugator. 

•  Phase-only  conjugator  (reflected  profile  is  flat-top). 

•  An  ordinary  mirror. 

By  comparison  of  the  beam  quality  with  the  conjugator  operating 
in  the  different  modes,  the  effect  of  the  conjugator  can  be 
clearly  demonstrated.  Simulations  of  a  phase-only  conjugator 
were  found  to  provide  useful  insight  into  the  RA/SBS  system  and 
were  a  good  complement  to  the  simulation  of  a  mirror  that 
reflects  the  correct  amplitude  but  the  wrong  phase.  This  condi¬ 
tion  could  be  representative  of  conjugation  by  stimulated  scat¬ 
tering  under  certain  conditions.23 

Nonlinear  propagation  is  accomplished  by  slicing  the  medium 
into  a  number  of  subregions.  The  nonlinear  interaction  between 
the  beams  is  calculated  at  the  surface  of  each  subregion  and  its 
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effects  are  incorporated  into  the  beams  prior  to  linear  propaga¬ 
tion  across  the  subregion.  The  accuracy  of  this  process  depends 
on  taking  an  adequate  number  of  subregions;  we  typically  used  60. 
The  nonlinear  interaction  takes  the  form 
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where  u  represents  the  complex  field,  k  the  propagation  constant, 
subscripts  p  and  s  the  pump  and  Stokes  waves,  and  7  =  70(l  +  i/3) 
the  complex  loss.  The  dispersion  factor  /?  is  included  to  model 
nonreciprocal  dispersion  effects,  such  as  when  the  seed  frequency 
is  off  line-center  of  the  Raman  gain  curve. 

As  discussed  above,  a  convergence  parameter  a  is  needed  to 
account  for  diffractive  size  changes  of  the  beams.  A  value  of  2 
was  found  necessary  to  give  adequate  beam  containment  within  the 
Raman  cell.  This  convergence  parameter,  however,  is  not  suffi¬ 
ciently  large  to  model  the  far-field,  as  is  demonstrated  on  the 
right  side  of  Figure  7.  As  seen  in  the  figure,  increasing  the 
convergence  parameter  to  a  =  8  gives  suitable  far-field  patterns. 
This  value  of  8,  however,  is  not  suitable  for  accurate  represen¬ 
tation  of  the  field  within  the  amplifier.  Incorporation  into  the 
code  of  a  different  convergence  transformation,  for  the  amplifier 
and  the  far  field  was  not  possible  within  the  time  span  of  this 
program.  As  a  result,  we  have  to  utilize  beam  quality  measure¬ 
ments  that  do  not  require  knowledge  of  the  full  extent  of  the 
beam  in  the  far-field.  The  Strehl  ratio  is  one  such  quantity;  it 
only  requires  an  accurate  calculation  of  the  peak  field  in  the 
f  ar-f ield . 

The  propagation  code  gives  the  intensity  and  phase  of  the 
seed  and  pump  at  any  point.  The  magnitude  or  phase  of  the  com¬ 
plex  mesh  point  values  are  output  as  a  grayscale  or  translated 
into  pseudo-color  plots.  Several  beam  quality  parameters  are 
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Figure  7.  Focal  plane  pattern  with  uncompensated 
aberrations  for  several  values  of  the 
convergence  parameter  a.  In  this  example, 
the  seed  wave  is  a  0.5-cm  flat-top  beam 
aberrated  by  a  MS  =  1.1  radian  phase 
aberrator,  without  Raman  amplification  or 
con j ugat ion . 

derived  from  these  mesh  point  values.  Photon  conversion  effi¬ 
ciency  is  calculated  by  integrating  the  intensity  of  the  seed  at 
the  end  of  the  Raman  cell  and  comparing  it  to  the  integrated 
input  pump  intensity.  The  near-field  beam  quality  is  evaluated 
using  the  phase-conjugate  fidelity: 

| /dxdy  u | uf | 2 

- (  (11) 

/ dxdy | u ;  | 2  / dxdy | u  f  |  2 

where  u  •,  and  Uf  are  the  initial  (prior  to  the  phase  aberration 
plate)  and  final  (after  conjugation,  gain,  and  passing  the  aber¬ 
ration  plate  a  second  time)  seed  (Stokes)  waves,  respectively. 
The  Strehl  ratio  is  calculated  from  the  peak  far-field  intensity 
pattern  and  is  equal  to  the  ratio  of  the  observed  peak  intensity 
and  the  peak  intensity  obtained  in  the  absence  of  aberrations 
(normalized  to  the  beam  energy) . 


A  comparison  of  the  conjugation  fidelity  with  the  Strehl 
ratio  for  an  assortment  of  geometries  is  shown  in  Figure  8.  As 
seen  in  the  figure,  the  two  measures  of  beam  quality  are  nearly 
identical.  In  no  case  did  we  find  that  the  two  measures  gave 
dramatically  different  results;  for  this  reason,  we  only  quote 
Strehl  ratios  in  this  section.  For  a  number  of  cases,  we  also 
calculated  the  far-field  beam  size  using  the  energy- in-a-bucket 
calculation.  A  representative  example  is  shown  in  Figure  9. 

Using  the  radius  that  contains  a  fixed  percentage  (84%  for  flat 
top  and  86%  for  Gaussian)  divided  by  the  diffraction-limited 
radius  gives  the  XDL  ( t imes-di f f ract ion- 1 imi ted)  of  the  beam. 
These  measurements  were  also  we  1 1 -cor re  1 ated  with  the  Strehl 
ratio.  That  is,  the  Strehl  ratio  is  the  inverse  of  the  XDL 
squared.  Once  again,  we  quote  only  the  Strehl  ratios  in  this 
section . 

As  described  above,  our  nonlinear  propagation  code  includes 
many  effects.  The  majority  of  these  effects  occur  because  of  the 
coupling  of  some  nonlinear  phenomena  (e.g.,  saturation)  with 
diffraction  in  the  presence  of  aberrations.  Diffractive  effects 
are  modeled  within  the  RA  as  well  as  outside  it.  It  follows  that 
a  great  number  of  effects  are  automatically  included  in  the  code. 
The  inclusion  of  Raman  gain  and  conservation  of  energy  allows  the 
observation  of  pump  replication  and  gain  saturation.  Nonlinear 
dispersion  was  included  by  allowing  a  complex  gain  coefficient. 
This  accounts  for  cases  where  the  SBS  frequency  shift  is  not 
fully  compensated.  The  inclusion  of  nonlinear  dispersion  causes 
large-  and  small-scale  inhomogeneities  in  the  phase  of  the  seed 
beam.  Large-scale  inhomogeneity  results  in  focusing  (or  defocus- 
ing)  of  the  seed  beam,  while  small-scale  inhomogeneities  result 
in  beam  breakup.  The  effects  of  RIVED  were  also  included  as 
discussed  below.  By  suitable  choice  of  the  degree  of  aberration, 
beam  geometry,  gain,  etc.,  we  were  able  to  observe  the  ability  of 
the  RA/SBS  system  to  compensate  for  the  reciprocal  phase  errors 
in  the  presence  of  the  nonreciprocal  phase  errors  discussed 
above . 
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Figure  8. 


Near-field  conjugation  fidelity  plotted 
as  a  function  of  the  far-field  Strehl  ratio  for 
an  assortment  of  geometries.  The  linear 
relationship  between  the  two  indicates  the 
correlation  of  these  two  measurements. 
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Several  additional  effects  could  have  been  incorporated  in 
the  code  if  time  had  permitted.  These  effects  include:  time- 
dependent  pulses,  very  high  Fresnel  number  geometries,  higher 
order  Stokes  waves,  and  broad  laser  linewidth  or  multiple  longi¬ 
tudinal  modes.  These  effects  were  not  included  because  they 
required  a  considerable  increase  in  the  execution  time,  and  the 
expense  was  not  judged  worthy.  Currently,  each  data  point 
requires  an  execution  time  of  about  9  min  on  the  DEC  VAX  8650,  or 
about  10  s  on  the  CRAY.  The  inclusion  of  time-dependent 
phenomena  would  require  us  to  store  the  values  in  each  subregion 
of  the  RA ;  this  amounts  to  a  large  increase  in  memory  require¬ 
ments  (~60X) .  High  Fresnel  number  geometries  would  require  a 
much  larger  mesh  grid,  smaller  subregion  thickness,  and  increased 
complexity  of  the  propagation  transformations.  Higher  order 
Stokes  waves  would  increase  the  complexity  in  proportion  to  the 
number  of  orders,  and  were  judged  to  be  not  important  in  the 
cases  we  studied.  Finally,  broad  linewidths  or  multiple  longi¬ 
tudinal  modes  would  also  increase  the  complexity  in  proportion  to 
the  number  of  modes,  and  were  judged  not  worthy  of  the  additional 
cost . 

3.2  GOALS  OF  THE  MODELING  PROGRAM/MODELING  TECHNIQUE 

The  goal  of  the  modeling  effort  was  to  determine  which  physi¬ 
cal  processes  limit  compensation  of  RAs  and  under  what  conditions 
they  do  so.  If  all  of  the  aberrations  were  reciprocal  and  if  the 
conjugation  fidelity  of  the  SBS  process  were  sufficiently  good, 
then  the  wavefront  of  the  amplified  field  should  be  very  nearly 
the  phase  conjugate  of  the  input  reference.  Departures  from  this 
ideal  performance  are  caused  by  nonreciprocal  aberrations  and 
imperfect  performance  of  the  SBS  conjugator.  The  purpose  of  this 
modeling  effort  was  to  quantify  the  amount  of  distortion  caused 
by  nonreciprocal  aberrations  for  uniform  as  well  as  distorted 
pump  and  seed  waves . 


We  looked  for  processes  that  deteriorate  the  beam  quality. 
These  processes  arise  from  imperfect  phase  conjugation,  the 
effects  of  nonlinear  gain,  diffraction,  and  the  mutual  inter¬ 
action  of  these  phenomena.  Deleterious  effects  from  the  conjuga- 
tor  can  arise  from  poor  replication  of  either  the  amplitude  or 
phase  of  the  input  seed.  The  effects  of  nonuniform  gain,  caused 
by  distortions  in  the  pump  and  saturation,  also  give  rise  to  a 
distorted  output.  Because  the  gain  in  the  RA  is  distributed,  the 
effect  of  nonuniform  gain  depends  on  the  scale  size  of  the  non- 
linearities,  relative  to  the  amount  of  diffraction  (pump  and/or 
the  seed) .  When  the  scale  size  of  the  gain  nonuniformity  is 
very  large  (uniform  gain),  the  RA/SBS  system  works  very  well.  As 
the  scale  size  decreases,  the  effect  of  nonuniform  gain  can 
dramatically  affect  the  fidelity  of  the  output.  In  the  limit  of 
very  short  scale  size,  the  fidelity  improves  due  to  averaging 
over  gain  nonuniformities.  It  was  found  that  the  conjugator  can 
compensate  for  reciprocal  phase  errors  in  almost  all  cases. 

The  advantage  of  utilizing  a  computer  model  to  investigate 
the  behavior  of  the  ampl i f i er /con j ugator  system  is  that  the 
various  nonreciprocal  aberrations  can  be  independently  intro¬ 
duced.  In  this  fashion,  the  deleterious  effects  of  each  can  be 
independently  evaluated.  It  was  not  possible,  nor  desirable,  to 
evaluate  all  combinations  of  effects  within  this  program. 

A  perfect  conjugator  can  compensate  for  phase  aberrations 
that  are  the  same  on  both  the  forward  and  backward  passes  (recip¬ 
rocal)  through  the  RA .  The  conjugator  can  not  be  expected  to 
compensate  for  aberrations  that  are  not  the  same  (nonreciprocal). 
Our  concern  is  that  nonlinear  gain  and  nonreciprocal  aberrations 
could  impair  the  ability  of  the  conjugator  to  compensate  for 
reciprocal  aberrations.  The  modeling  part  of  this  study  has  been 
directed  towards  either  strengthening  or  rejecting  this  concern. 

We  chose  to  quantify  these  results  by  asking  the  following 
questions:  Under  what  conditions  are  strong  nonlinear  effects 

expected  to  arise0  Under  these  conditions,  How  does  the  beam 
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quality  of  a  Stokes  wave  derived  from  an  aberrated  seed  and 
reflected  from  a  conjugate  mirror  compare  to  the  beam  quality  of 
a  Stokes  wave  derived  from  a  clean  seed  and  reflected  from  an 
ordinary  mirror?  The  beam  quality,  calculated  from  the  far-field 
Strehl  ratio  or  the  near-field  conjugation  f i .  lity,  will  be  poor 
in  either  case  if  the  nonreciprocal  errors  are  strong.  However, 
a  comparison  of  the  aberrated-con j ugated  wave  with  the  unaber¬ 
rated  wave  highlights  the  way  in  which  the  compensation  for 
reciprocal  errors  is  jeopardized  by  the  Raman  amplification 
process . 

As  described  above,  the  nonlinear  propagation  code  allows 
variation  of  a  number  of  parameters,  including: 

•  Pump  wavelength 

•  Stokes  wavelength 

•  Beam  diameters 

•  Raman  cell  length 

•  Small-signal  gain  parameter  (gIL) 

•  Initial  seed-to-pump  amplitude  ratio  (seed  intensity) 

•  RMS  phase  error  of  aberration  plates 

•  Outer  scale  length  of  aberration  plates 

•  Choice  of  amplitude  profile  for  pump  and  seed  beams 

•  Choice  of  conjugator: 

perfect  phase  and  amplitude  conjugator 
phase-only  conjugator 
mirror  reflector. 

Because  of  the  large  number  of  parameters,  it  was  impossible 
to  investigate  all  possible  values  within  the  scope  of  this  pro¬ 
gram.  In  fact,  not  all  combinations  of  parameters  are  interest¬ 
ing.  In  order  to  make  the  modeling  effort  tractable,  we 
attempted  to  limit  our  investigation  to  combinations  of  para¬ 
meters  that  provide  useful  insight  into  the  operation  of  the 
RA/SBS  system.  We  chose  four  standard  beam  profiles  for  modeling 


(1)  3 . O-cm-diameter  flat-top  beam. 

(2)  0 . 5-cm-diameter  flat-top  beam. 

(3)  Gaussian  beam,  truncated  at  0 . 75-cm-diameter . 

(4)  Gaussian  beam,  truncated  at  0 . 35-cm-diameter . 

The  1/e  radius  of  the  Gaussian  intensity  profile  is  set  to  2/3  of 
the  truncation  radius.  For  beam  profiles  (1)  and  (2),  the 
injected  seed  intensity  was  10‘4  of  the  pump  intensity;  for  beam 
profile  (3)  it  was  10* 2  of  the  pump  intensity,  and  for  profile 
(4)  it  was  0.16  of  the  pump  intensity.  Raman  cell  length  was 
held  at  100  cm.  Beam  profile  (1)  corresponds  to  a  Fresnel  number 
(a2/LX)  of  ~350  and  represents  the  largest  Fresnel  number  that 
can  be  comfortably  handled  by  the  code.  High  Fresnel  numbers  are 
commonly  encountered  in  high  energy  applications  where  optical 
damage  thresholds  place  a  lower  limit  on  beam  size.  Larger 
Fresnel  numbers  require  a  large  increase  in  calculation  time,  as 
discussed  above.  The  Fresnel  number  of  beam  profile  (2)  is  ~9 . 1 
and  is  the  lowest  Fresnel  number  discussed  in  this  section.  Beam 
profile  (3)  has  an  intermediate  Fresnel  number  (~21).  In  addi¬ 
tion,  this  beam  is  a  truncated  Gaussian,  resulting  in  a  decreased 
far-field  diffraction  angle  and  increased  effects  from  nonuniform 
intensity  profiles.  Finally,  profile  (4),  with  a  Fresnel  number 
of  ~4.5,  closely  corresponds  to  our  experiments  using  a 
frequency-doubled  single-mode  Nd:YAG  laser. 

After  choosing  the  beam  geometries  discussed  above,  it  was 
also  necessary  to  determine  how  to  incorporate  the  nonlinear 
effects.  The  difficulty  here  was  to  discriminate  between  the 
individual  effects.  To  gain  insight  into  the  nonlinear  effects, 
we  utilized  the  following  sequence: 

(a)  Seed  aberration,  diffraction  and  conjugation. 

(b)  Add  the  effects  of  Raman  gain  to  case  (a) . 

(c)  Add  the  effects  of  pump  aberration  to  case  (b) . 

(d)  Add  the  effects  of  RIVED  to  case  (c)  . 

(e)  Add  the  effects  of  nonlinear  dispersion  to  case  (b) . 


We  found  this  sequence  provided  insight  into  the  ability  of  the 
conjugator  to  compensate  for  reciprocal  phase  errors  under 
increasingly  difficult  situations.  The  results  of  this  investi¬ 
gation  are  discussed  in  detail  in  the  following  subsection. 

3.3  RESULTS 

The  results  of  this  analysis  are  summarized  below  with  refer¬ 
ence  to  Figures  10  to  17.  Simulations  of  our  experimental  pro¬ 
gram  are  reviewed  in  Section  5.  Our  results  generally  consist  of 
the  Strehl  ratio  shown  by  curves  representing  the  three  experi¬ 
mental  conditions:  first,  an  unaberrated  seed  beam  and  a  simple 
reflection  rather  than  phase  conjugation;  second,  a  highly  aber¬ 
rated  seed  and  a  perfect  phase  and  amplitude  conjugator;  third,  a 
highly  aberrated  seed  and  a  phase  only  conjugator  (the  conjugate 
amplitude  is  a  flat-top) .  The  Strehl  ratio  in  these  three  condi¬ 
tions  is  plotted  as  we  add  in  nonlinear  and  nonreciprocal 
effects.  By  comparing  the  unaberrated  results  with  the 
aberrated-con j ugated  results,  we  will  be  able  to  identify  those 
conditions  under  which  the  conjugator  is  unable  to  compensate  for 
reciprocal  phase  errors  in  the  presence  of  nonreciprocal  effects. 
Comparison  of  the  perfect  conjugation  results  with  the  phase-only 
conjugation  results  indicates  the  importance  of  high  fidelity 
replication  of  the  intensity  profile  during  conjugation.  That 
is,  some  of  the  nonreciprocal  effects,  such  as  gain  saturation, 
alter  the  intensity  profile  in  a  manner  equivalent  to  the  phase- 
only  conjugation. 

To  begin  this  analysis,  we  first  examine  the  effects  of  seed 
aberration  -  that  is,  reciprocal  phase  errors.  In  Figure  10, 
three  sets  of  data  are  presented;  these  correspond  to  intensity 
profiles  (1),  (2),  and  (3)  discussed  above.  In  the  upper  left 

quadrant  of  the  figure,  the  high  Fresnel  number  (350),  3.0  cm 
flat-top  beam  results  are  given.  The  upper  right  quadrant  con¬ 
tains  the  results  for  the  low  Fresnel  number  (9.1)  flat-top  bea"i. 
The  lower  left  hand  quadrant  gives  the  results  for  the  medium 


STREHL  RATIO  STREHL  RATIO 


STREHL  RATIO  STREHL  RATIO 


0.6  - 

0.4  - 

0.2  - 

0.0  - 


0.75  cm  GAUSSIAN 


Figure  11.  The  effect  of  Rama 
ratio  of  the  seed 
top  curve  is  for  a 
the  central  curve 
conjugated  geometr 
the  aberrated-phas 


ft  X  h-  3  o'  3 


STREHL  RATIO 


0.75  cm  GAUSSIAN 


INJECTION  RATIO  =  10'1 


INJECTION  RATIO  =  10"2 


INJECTION  RATIO  =  10' 


INJECTION  RATIO  =  10' 


A0S  =8.2  rad 


0.0  0.1 


Figure  12. 


0.2  0.3  0.4  0.5  0.6  0  7  0.8  0  9  1 

PHOTON  CONVERSION  EFFICIENCY 


The  effect  of  varying  the  injection 
ratio  on  the  Strehl  ratio  of  the  seed  beam 
as  a  function  of  the  photon  conversion 
efficiency.  In  this  example,  0.75-cm- 
diameter  truncated  Gaussian  beams  are  used 
with  an  RMS  seed  phase  aberration  of  8.2  rad 
and  no  pump  aberration. 


WWW 


aberrated  with  A0RMS  = 
1-signal  gain  paramete 


■  »  -  r  -  ir-  ir^  vmvrv.  r*.  ST»  VKVWV 


3.0  cm  FLAT-TOP 


0.5  cm  FLAT-TOP 


0.75  cm  GAUSSIAN 


0.5  cm  FLAT-TOP 


A0S=  8.2 


gIL  =  8.0 


RMS  PUMP  PHASE  ERROR,  rad 


RMS  PUMP  PHASE  ERROR,  rad 


Figure  15. 


The  effect  of  pump  aberration  on  the  Strehl 
ratio  of  the  seed  beam.  In  each  quadrant,  the 
top  curve  is  for  an  unaberrated  beam  geometry, 
the  central  curve  the  aberrated-con j ugated 
geometry,  and  the  lowest  curve  the  aberrated- 
phase-only-conj ugated  geometry. 
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Fresnel  number  (21)  Gaussian  beam.  In  each  quadrant,  we  plot  the 
Strehl  ratio  of  the  seed  beam  as  a  function  of  the  seed  aberra¬ 
tion.  The  amount  of  seed  aberration  is  represented  by  the  RMS 
phase  error  of  the  phase  aberration  plate.  This  RMS  value  is 
found  by  direct  evaluation  of  the  phase  screen,  because  different 
initialization  values  for  the  random  number  generator  yield 
slightly  different  RMS  phase  errors.  RMS  phase  errors  up  to 
10  rad  are  shown  in  the  figure.  In  RAs  where  the  gas  flow  has 
been  specifically  designed  to  minimize  thermal  gradations,  we 
expect  aberrations  of  a  few  radians  (see  Appendix  B) .  It  should 
be  remembered  in  the  following  discussion  that  we  are  considering 
optical  beams  aberrated  up  to  2  to  3  times  what  would  be  expected 
in  a  well-designed  system.  In  the  absence  of  phase  conjugation, 
the  aberration  of  the  seed  beam  is  actually  twice  that  of  the 
aberrating  plate  because  it  is  double  passed  in  the  system. 

The  top  curve  in  each  quadrant  gives  the  Strehl  ratio  of  the 
output  beam  when  a  perfect  phase  and  amplitude  phase  conjugator 
is  used.  As  expected,  there  is  no  reduction  of  Strehl  ratio  for 
RMS  phase  aberrations  as  large  as  8  rad.  The  figures  show  a 
Strehl  ratio  of  0.985  for  the  flat-top  beams;  this  represents  the 
precision  with  which  our  code  can  numerically  diffract  beams. 

With  the  truncated  Gaussian  beam,  the  results  were  normalized  to 
the  unaberrated  case.  As  a  result,  the  Strehl  is  1.0  for  the 
perfect  conjugator.  The  Strehl  of  the  conjugated  beam  remains  at 
1.0  for  seed  aberrations  as  large  as  11  rad,  demonstrating  the 
ability  of  our  diffraction  code  to  accurately  simulate  diffrac¬ 
tion  of  highly  aberrated  beams.  The  central  line  in  each  quad¬ 
rant  represents  the  case  where  a  phase-only  conjugator  is  used. 

As  expected,  the  Strehl  ratio  falls  off  with  increasing  aberra¬ 
tion.  For  small  aberrations,  the  fidelity  of  the  flat-top  beams 
is  very  high,  because  propagation  across  the  amplifier  cell  has 
not  turned  the  small  phase  aberration  into  a  significant  ampli¬ 
tude  aberration.  Therefore,  the  phase-only  (flat-top  amplitude) 
conjugation  gives  a  reasonably  good  conjugate  of  the  input  beam. 


It  is  evident  that  the  beam  fidelity  falls  off  with  increas¬ 
ing  aberration  because  diffractive  effects  induce  intensity 
variations  on  the  beam  before  it  reaches  the  conjugator.  It  is 
interesting  to  note  that  the  Strehl  ratio  falls  off  much  more 
rapidly  for  the  low  Fresnel  number  beams,  where  diffractive 
effects  are  much  stronger.  For  small  aberrations,  the  Strehl 
ratio  of  the  Gaussian  beam  does  not  asymptotically  approach  the 
perfect  conjugator  case,  because  we  assume  a  flat-top  intensity 
profile  for  our  phase-only  conjugator  even  for  the  Gaussian  beam 
case.  Finally,  the  lower  curve  in  each  quadrant  gives  the 
results  when  an  ordinary  mirror  is  used.  As  can  be  seen,  the 
Strehl  of  the  unconjugated  system  falls  off  very  rapidly,  reduc¬ 
ing  to  ~0 . i  with  RMS  phase  errors  of  ~1.3  rad.  In  the  descrip¬ 
tion  that  loilows,  we  find  interesting  results  for  RA/SBS  systems 
where  the  RMS  phase  errors  are  ~8  rad.  Here,  the  uncorrected 
(i.e.,  without  conjugation)  Strehl  ratio  is  <<0.001;  the  actual 
value  in  this  case  is  subject  to  significant  numerical  errors  and 
is  too  small  to  be  of  significance.  For  this  reason,  we  do  not 
attempt  to  compare  the  conjugated  system  with  the  unconjugated 
system;  rather,  the  conjugated  system  with  seed  phase-aberrations 
is  compared  to  an  unconjugated  system  without  seed  phase-aberra¬ 
tions.  It  should  be  remembered  throughout  this  discussion  that, 
even  when  the  Strehl  ratio  of  the  conjugated  system  drops  by  an 
order  of  magnitude,  we  are  modeling  beams  with  such  large  aberra¬ 
tions  that,  without  conjugation,  the  Strehl  ratio  would  be  ~103 
lower . 

Figure  11  shows  the  effect  of  adding  gain  to  the  RA .  To 
demonstrate  the  effect  of  gain  in  the  amplifier,  we  have  chosen 
four  representative  aberration  conditions,  indicated  by  the 
arrows  on  Figure  10.  Three  curves  are  shown  in  each  quadrant  of 
the  figure.  The  upper  curve  gives  the  output  Strehl  ratio  of  an 
unaberrated  seed  beam.  The  central  line  shows  the  results  for  a 
aberrated-con j ugated  beam.  The  lower  line  represents  an 
aberrated -phase -on  1 y-con j ugated  beam.  In  all  cases  the  pump  beam 


was  unaberrated.  With  very  low  Raman  gain,  the  perfect  conjuga- 
tor  is  able  to  correct  for  all  seed  aberrations,  while  the  phase- 
only  conjugator  is  somewhat  less  effective.  As  the  gain  is 
increased,  saturation  and  pump  nonuniformity  begin  to  play  an 
important  role,  altering  the  intensity  pattern  of  the  seed  beam 
on  its  second  pass  through  the  amplifier.  This,  together  with 
diffraction,  yields  an  inaccurate  intensity/phase  profile  when 
the  beam  retraverses  the  seed-phase  aberrating  plate,  and  leads 
to  a  reduced  Strehl  ratio.  When  the  gain  becomes  extremely 
high,  the  pump  is  strongly  depleted  soon  after  it  enters  the 
amplifier  (before  diffraction  has  a  significant  effect).  This 
results  in  a  seed  amplitude  profile  that  nearly  matches  the  pump 
amplitude  profile  at  the  amplifier  input.  Because  the  pump  is 
unaberrated  in  this  example,  the  seed  profile  becomes  uniform 
(flat-top),  closely  matching  the  case  of  a  phase-only  conjugator. 
As  can  be  seen  in  the  figure,  the  Strehl  ratio  of  the  aberrated- 
conjugated  beam  matches  the  unaberrated  beam  for  low  gain,  and 
approaches  the  Strehl  ratio  of  the  phase-only  conjugated  beam  at 
high  gains. 

The  results  for  the  low  Fresnel  number  flat-top  beam  are 
shown  on  the  right  side  of  Figure  11  In  the  upper  right  quad¬ 
rant  are  the  results  for  low  seed  aberration  (#RMS  =1-4  rad). 

We  see  that,  for  typical  values  of  small-signal  gain  (gIL  <10) 
the  phase  conjugator  is  capable  of  correcting  for  essentially  all 
the  reciprocal  phase  errors  (induced  by  the  aberration  plate). 
This  is  evident  because  the  Strehl  ratio  of  the  aberrated-con- 
j ugated  beam  is  indistinguishable  from  the  unaberrated  beam. 

Both  Strehl  ratios  are  reduced  because  the  diffraction  rings  on 
the  pump  lead  to  nonuniform  gain,  altering  the  intensity  profile 
of  the  seed  beam.  At  high  gains  (gIL  ~14) ,  pump  depletion  begins 
to  play  an  important  role.  The  Strehl  ratio  of  the  unaberrated- 
unconjugated  beam  continues  to  decrease,  approaching  a  value  of 
~0 . 5 .  The  Strehl  ratio  of  the  aberrated-conj ugated  beam 
decreases  only  slightly  farther  and  becomes  equaj  to  the  phase- 
only-conjugation  Strehl  at  very  high  gains  (saturation). 


A  more  dramatic  example  is  shown  in  the  lower  right  quadrant 
where  the  seed  phase  aberration  has  been  increased  ~6X 
($RMS  =  8.2  rad) .  The  general  trends  are  the  same  as  discussed 
above;  however,  the  effect  of  aberration  is  much  stronger.  The 
difference  is  the  intensity  aberration  of  the  seed  (due  to  dif¬ 
fraction  of  larger  seed  aberrations)  is  much  larger,  and  the 
effects  of  poor  amplitude  preservation  (gain  saturation  or  phase- 
only  conjugation)  are  more  evident. 

As  discussed  above,  the  decrease  of  the  Strehl  ratio  under 
conditions  of  high  gain  is  due  largely  to  intensity  aberrations 
that  result  from  diffraction  of  the  large  phase  aberrations  pro¬ 
duced  by  the  aberration  plate.  It  follows  that  a  higher  Fresnel 
number  geometry  should  yield  better  results.  The  high  Fresnel 
number  beam  geometry  is  given  in  the  upper  left  quadrant  of 
Figure  11.  As  in  the  lower  right,  the  phase  aberration 
0rms  =  8.2  rad.  Here,  the  Strehl  ratio  remains  largely  indepen¬ 
dent  of  the  Raman  gain,  because  the  large  Fresnel  number  geometry 
minimizes  diffraction  effects  and  the  seed  beam  intensity  pro¬ 
file  remains  uniform. 

The  Gaussian  beam  result  is  shown  in  the  lower  left  quadrant 
of  Figure  11.  The  results  are  qualitatively  different  from  those 
obtained  with  flat-top  beams:  at  low  to  moderate  gains  the 

Strehl  ratio  decreases,  but  at  high  gains  the  ratio  increases. 
Investigation  has  shown  that  this  initial  decrease  is  due  to  the 
nonuniform  gain  resulting  from  the  Gaussian  profile  of  the  pump, 
which  serves  to  dramatically  reduce  the  near-field  beam  diameter 
and  hence  the  Strehl  ratio.  At  high  gains,  saturation  becomes 
important.  Because  the  gain  is  significantly  higher  in  the  cen¬ 
ter  of  the  pump,  the  seed  beam  depletes  the  central  portion  of 
the  pump,  giving  relatively  higher  gain  to  the  outer  regions.  As 
the  seed  depletes  the  entire  pump  beam,  it  begins  to  replicate 
the  pump  intensity  profile.  If  it  were  not  for  diffraction,  the 
two  intensity  profiles  would  be  identical  and  the  Strehl  ratio 
would  regain  a  value  very  near  1 . 


To  emphasize  this  point,  we  have  plotted  the  aberrated-con- 
jugated  curve  for  the  Gaussian  beam  along  with  similar  curves 
calculated  for  other  injection  ratios  (seed-to-pump  intensity 
ratio) .  As  the  injection  ratio  changes,  so  does  the  amount  of 
gain  required  for  the  seed  beam  to  deplete  the  pump.  To  compen¬ 
sate  for  this,  we  plot  the  data  in  Figure  12  as  a  function  of 
photon  conversion  efficiency  (rj)  rather  than  the  small-signal 
gain  parameter  (gIL) .  The  aberrated-con j ugated  curve  of  the 
lower  left  quadrant  of  Figure  11  is  shown  by  the  curve  labeled 
"injection  ratio  =  10'2";  that  is,  an  injection  ratio  of  1%.  It 
is  clear  from  the  figure  that,  as  the  injection  ratio  is  reduced, 
the  Strehl -reduci ng-ef f ect  of  nonuniform  gain  is  enhanced.  It  is 
interesting  that  the  Strehl  ratio  for  lower  injection  remains 
low,  even  at  high  photon  conversion  efficiency.  The  reason  is 
the  low  injection  ratio  requires  very  high  gain  and,  with  high 
gain,  aberration  of  the  pump  becomes  very  important.  For  the 
purposes  of  this  plot,  we  are  using  an  RMS  seed  phase  error  of 
8.2  rad.  When  the  unconjugated  seed  double  passed  this  aberra- 
tor,  its  Strehl  ratio  was  reduced  to  well  below  10'3.  In  fact, 
to  the  extent  that  single-passing  an  aberrator  with 
A0RMS  ~8 . 2  rad  is  equivalent  to  double-passing  a  A$RMS  ~4 . 1  rad 
aberrator,  we  expect  that  the  Strehl  ratio  of  the  pump  is  ~0.01 
(see  Figure  10)  after  passing  through  the  seed  aberration  plate. 
Therefore,  if  the  seed  beam  truly  matched  the  pump  (pump  replica¬ 
tion),  the  Strehl  ratio  would  be  ~0.01,  rather  than  the  observed 
values  >0.3.  It  follows  that,  even  in  the  worst  case  represented 
in  Figure  12,  the  conjugator  still  does  a  good  job  of  correcting 
the  reciprocal  phase  errors  . 

In  Figure  13,  we  have  plotted  the  photon  conversion  effi¬ 
ciency  as  a  function  of  the  small-signal  gain  parameter.  Shown 
in  the  figure  are  the  results  for  a  0 . 5-cm-di ameter  flat-top  beam 
(injection  ratio  10'4)  and  a  0.75-cm  truncated  Gaussian  beam 
(injection  ratios  10‘2  and  10‘4).  It  is  clear  from  the  figure 
that  it  is  much  more  difficult  to  deplete  the  Raman  gain  (high 


photon  conversion  efficiency)  for  Gaussian  beams  than  it  is  for 
flat-top  beams.  The  flat-top  beam  attains  good  photon  conversion 
efficiency  (>90%)  with  physically  realizable  small-signal  gain 
parameters  (gIL  ~15) .  This  flat-top  evaluation  is  in  good  agree¬ 
ment  with  the  formula  for  photon  conversion  efficiency  ( rj )  in  one 
dimension . 2  4 
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The  equivalent  situation  for  a  Gaussian  beam  is  shown  by  the 
lowest  curve  in  Figure  13.  At  a  small-signal  gain  parameter 
(gIL)  of  15,  the  photon  conversion  efficiency  of  the  Gaussian 
beam  is  roughly  half  that  of  the  flat-top  beam.  At  the  highest 
gains  shown  in  the  figure,  the  photon  conversion  efficiency  has 
still  attained  only  a  value  of  ~60%  —  the  result  of  the  low  gain 
on  the  outer  edges  of  the  pump.  If  the  injection  ratio  of  the 
seed  is  increased,  the  situation  improves.  The  central  curve  is 
plotted  for  a  Gaussian  beam  with  an  injection  ratio  100X  larger 
than  the  other  two  curves.  For  small  gains,  the  conversion  effi¬ 
ciency  is  hi  gher  ,  but  the  flat-top  beam  is  superior  for  high 
conversion  efficiency.  Analytic  results  for  small-signal  gain 
parameters  as  high  as  25  are  plotted  in  the  figure.  In  practice, 
however,  these  gains  are  probably  not  achievable.  Not  only  will 
amplified  spontaneous  emission  (ASE)  become  important,  but  also 
2nd  Stokes  will  be  generated  from  the  center  before  the  outer 


rang  is  s 


aturated .  The 


:ffects  were  not  included  in  the  code. 


Up  to  this  point,  we  have  discussed  the  effects  of  reciprocal 
phase  error  and  diffraction  in  the  RA  with  a  unaberrated  pump 
beam.  The  effects  of  depletion  and  saturation  of  the  pump  beam 
have  also  been  discussed.  When  considering  a  truncated  Gaussian 
pump  beam,  we  found  significant  effects  caused  by  nonuniform 
gain.  At  this  point,  we  wish  to  consider  the  effects  of  adding 
phase  and  intensity  aberrations  to  the  pump  beam,  which  results 
in  highly  nonuniform  gain  that  cannot  be  compensated  with  the 
phase  conjugate  mirror. 


We  created  a  dirty  pump  beam  with  both  phase  and  amplitude 
aberratior  by  placing  a  phase  aberrating  screen  in  the  path  of 
the  pump  learn  and  allowing  propagation  (diffraction)  to  create 
intensity  aberrations.  For  convenience,  this  propagation  dis¬ 
tance  was  chosen  to  be  equal  to  the  Raman  cell  length.  The 
intensity  aberrations  of  the  pump  beam  are  transmitted  to  the 
seed  beam  through  nonuniform  gain.  These  seed  intensity  aberra¬ 
tions  develop  into  phase  aberrations  during  propagation  across 
the  amplifier.  The  degree  to  which  the  seed  beam  replicates  the 
pump  depends  on  the  interplay  of  nonuniform  gain,  saturation,  and 
diffraction.  High  fidelity  replication  of  the  pump  occurs  only 
under  limited  conditions  and  was  not  observed  from  the  simula¬ 
tions  done  under  this  program.  An  example  of  this  process  is 
shown  in  Figure  14.  On  the  left  side  of  the  figure,  the  seed 
beam  (a  0.5-cm  unaberrated  flat-top)  is  shown  at  the  entrance  to 
the  amplifier  (after  conjugation).  The  upper  left  grayscale  plot 
is  the  intensity  profile,  while  the  lower  left  grayscale  plot  is 
the  phase  profile.  It  is  evident  that  the  seed  has  developed  a 
small  phase  ripple  as  a  result  of  diffraction,  but  there  is  no 
discernible  ripple  in  the  intensity.  The  pump  beam,  shown  by  the 
central  grayscale  plots,  has  been  aberrated  (A0 RMS  =  3.8  rad)  and 
allowed  to  propagate.  The  pump  has  large  amplitude  and  phase 
aberrations  (pump  Strehl  ~0.025) .  On  the  right  side  of  the 
figure,  the  intensity  and  phase  of  the  seed  are  shown  at  the 
amplifier  output  (gIL  =  8) .  Clearly,  the  seed  has  picked  up 
large  aberrations  from  the  pump.  The  seed  aberrations  do  not, 
however,  replicate  the  pump  aberrations  because:  (1)  the  pump  and 
seed  diffract  at  different  rates  because  of  their  different  wave¬ 
lengths;  (2)  the  small-signal  gain  parameter  gIL  =  8  causes  the 
high  intensity  portions  of  the  pump  to  saturate  at  the  beginning 
of  the  RA ;  consequently,  nonuniform  gain  and  diffraction  can  not 
cooperate  to  give  a  good  replication. 

For  comparison  with  prior  results,  we  have  considered  the 
addition  of  pump  aberrations  to  the  results  shown  in  Figure  11. 
Figure  15  shows  the  effect  of  aberrating  the  pump  beam  at  the 
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points  indicated  by  the  arrows  in  Figure  11.  The  beam  geometries 
for  each  quadrant  in  Figure  15  are  the  same  as  in  the  correspond¬ 
ing  quadrant  in  Figure  11.  The  general  conclusion  from  this 
figure  is  that  nonuniform  gain  caused  by  pump  intensity  aberra¬ 
tions  cannot  be  compensated  by  phase  conjugation.  This  does  not 
mean  that  the  conj  ugator  is  not  working.  The  conjugator  is  able 
to  compensate  for  essentially  all  the  reciprocal  phase  errors 
introduced  by  the  aberration  plate.  It  is  clear  from  Figure  15 
that  the  Strehl  ratio  for  all  three  curves  fall  off  together. 

That  is,  the  aber;  tei-con  j  ugated  system  behaves  just  as  well  as 
the  unaberrated  system.  It  follows  that  gain  nonuniformities 
must  be  minimized  by  "smoothing"  out  the  pump  intensity  aberra¬ 
tions.  This  can  be  accomplished  by  providing  more  uniform  pump 
beams  or  by  lowering  the  gain  and  increasing  aberration  so  that 
the  pump  beam  aberration  averages  out.  It  is  also  evident  from 
Figure  15  that  the  higher  Fresnel  number  geometries  (large  flat- 
top  and  Gaussian  beams)  have  higher  Strehl  ratios,  which  is  not 
surprising,  because  the  larger  beams  experience  less  diffraction 
and  hence  less  intensity  aberration. 

As  described  in  the  preceding  discussion,  we  found  no  situa¬ 
tions  under  which  the  perfect  phase  and  amplitude  conjugator  was 
unable  to  compensate  for  essentially  all  the  reciprocal  phase 
errors.  We  do,  however,  find  such  an  effect  when  we  consider 
RIVED,  which,  in  contrast  to  the  effects  described  above,  leads 
to  a  direct  phase  error  rather  than  an  intensity  error. 

The  Raman  amplification  process  transfers  molecules  from  the 
ground  state  to  a  vibrationally  excited  state.  While  in  this 
excited  state,  the  molecules  have  a  different  polarizability, 
which  leads  to  a  change  in  the  index  of  refraction.  Clearly, 
this  can  only  be  a  strong  effect  when  significant  numbers  of 
molecules  are  transferred  to  the  excited  state,  and  when  they  do 
not  have  time  to  thermalize  and/or  diffuse  out  of  the  regions  of 
high  intensity.  If  all  the  molecules  are  excited,  the  polari¬ 
zability  of  the  medium  increases  by  20%  (hydrogen) .  To  estimate 
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the  change  in  refractive  index,  we  assume  an  optical  pulse  that 
is  much  shorter  than  the  phonon  relaxation  time.  It  follows  that 
every  photon  transferred  from  the  pump  to  the  Stokes  wave  yields 
an  excited  molecule.  If  we  also  assume  that  the  medium  does  not 
saturate  (i.e.,  we  don’t  use  up  the  available  molecules),  we  find 
that 
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(13) 


where  A<I>  is  the  phase  error  accumulated  by  the  end  of  the  pulse 
at  a  given  location,  P  the  gas  pressure,  (n-l)0  ~135  x 
10"6  amagat'1,  (Jp/A)  the  initial  peak  pump  intensity,  N  the  gas 
density,  h  is  Planck’s  constant,  c  the  speed  of  light,  and  us  the 
Stokes  field  normalized  (as  in  the  code)  to  the  initial  peak  pump 
field. 

The  result  of  this  simulation  is  shown  in  Figure  16.  As  in 
Figures  11  and  15,  there  are  four  quadrants,  each  representing  a 
particular  beam  geometry.  The  upper  and  lower  right  quadrants 
are  for  the  low  Fresnel  number  flat-top  beam,  with  the  low  aber¬ 
ration  at  the  top  and  high  aberration  at  the  bottom.  The  upper 
left  is  for  a  high  Fresnel  number  flat-top  beam,  and  in  the  lower 
left,  a  moderate  Fresnel  number  truncated  Gaussian.  Each  quad¬ 
rant  contains  four  curves.  These  curves  represent  the  unaber¬ 
rated  and  the  aberrated-con j ugated  beams  with,  and  without,  pump 
aberration;  each  is  plotted  as  a  function  of  pump  energy  density. 
The  unaberrated  results  are  indicated  by  the  solid  circles,  while 
the  aberrated-con j ugated  results  are  shown  by  the  solid  squares. 

As  expected,  the  Strehl  ratio  is  reduced  by  the  effects  of 
RIVED.  It  is  interesting  that,  when  the  pump  is  unaberrated 
(A0p  =  0.0  rad),  the  Strehl  ratio  of  the  aberrated-con j ugated 
beam  can  be  significantly  lower  than  that  of  the  unaberrated- 
unconjugated  beam.  In  contrast,  the  two  Strehl  ratios  are  nearly 
equivalent  for  aberrated  pump  beams. 
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For  unaberrated  pump  beams,  we  believe  the  reduction  in 
Strehl  ratio  is  due  to  the  effect  of  seed  aberration.  That  is, 
the  seed  beam  aberration  causes  the  seed  intensity  near  the  con- 
j ugator  to  be  highly  nonuniform.  The  number  of  molecules  trans¬ 
ferred  to  the  excited  vibrational  state,  giving  rise  to  the 
RIVED,  is  then  proportional  to  the  seed  intensity  (uniform  pump 
intensity),  and,  thus,  is  spatially  highly  nonuniform.  It 
follows  that  the  aberrated  seed  beam  experiences  a  much  larger 
nonreciprocal  phase  error.  For  the  Gaussian  beam,  the  effect  is 
even  more  pronounced  —  the  result  of  a  lensing  action  caused  by 
to  the  higher  pump  intensity,  and  therefore  higher  photon  conver¬ 
sion  efficiency,  near  the  center  of  the  beam. 

For  the  aberrated  pump  beam,  spatial  nonuniformity  of  the 
excited  vibrational  state  density  results  from  pump  intensity 
aberrations  as  well  as  seed  aberrations.  In  this  case,  the 
effect  from  the  seed  aberration  is  masked  by  the  pump  aberration. 

This  RIVED  effect  is  applicable  only  for  short  optical 
pulses  -  so  short  that  the  seed  beam  does  not  overlap  itself 
completely  on  the  first  and  second  passes  through  the  amplifier. 
If  it  did,  the  phase  error,  due  to  RIVED,  would  be  sensed  on  the 
first  pass  and  compensated  by  the  conjugation  (i.e.,  the  effect 
i s  reciprocal ) . 

An  additional  effect,  giving  rise  to  a  direct  phase  aberra¬ 
tion  of  the  seed  beam,  is  due  to  nonlinear  dispersion,  which  in 
turn  is  a  result  of  frequency  mismatch  between  the  Stokes  seed 
beam  and  the  Raman  gain.  This  effect  can  result  from  either 
multimode  operation,  where  the  laser  linewidth  is  not  narrow 
compared  to  the  Raman  linewidth,  or  from  an  uncompensated  SBS 
frequency  shift.  In  contrast  to  RIVED,  this  is  a  true  nonrecip¬ 
rocal  phase  error,  because  the  Raman  process  gives  gain  only  in 
one  direction.  As  discussed  in  Section  2,  the  amount  of  non¬ 
linear  dispersion  can  be  evaluated  from  the  Kramers-Kronig  rela¬ 
tionship.  The  magnitude  of  the  dispersion  is  proportional  to  the 
small-signal  gain  of  the  amplifier: 
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gIL  =  -  g0IL  .  (14) 

1  +  4  (5w/T ) 2 

This  relationship  shows  that  the  maximum  dispersion  occurs  at  a 
frequency  offset  equal  to  one-half  the  f ul 1 -wi dth-hal f -maximum 
(Su  =  T/2)  of  the  Raman  gain  line;  at  this  point,  it  attains  a 
value  of  g0IL/2.  Nonlinear  dispersion  has  been  simulated  in  our 
code  by  using  a  complex  gain  parameter  gIL  =  <t>(l  -  i/c)  .  It  fol¬ 
lows  that  the  largest  realizable  value  of  K  is  0.5.  From 
Figure  17,  which  shows  the  results  of  modeling  this  effect,  the 
conclusion  can  be  drawn  that  this  value  will  not  cause  any 
significant  degradation  of  beam  quality. 

It  follows  that  nonlinear  dispersion  does  not  have  a  serious 
effect  on  the  RA/SBS  system.  For  illustration  purposes,  we  have 
modeled  frequency  offsets  as  large  as  K  =  10.  Although  this 
represents  a  condition  not  achievable  by  the  process  of  nonlinear 
dispersion,  there  could  be  other  nonlinear  processes  that  pro¬ 
duce  such  effects.  The  figure  shows  that  the  effect  is  more 
significant  for  distorted  pump  beams  and  can  have  significant 
effect  for  values  of  k  approaching  2.  This  implies  that  the 
nonlinear  induced  phase  shift  must  be  kept  to  less  than  a  wave 
within  the  Raman  gain  length  (i.e.,  K  <1). 

In  summary,  we  have  modeled  the  effects  of  pump  and  seed 
aberration  on  the  beam  fidelity  of  a  RA/SBS  system.  Our  simula¬ 
tion  includes  diffraction  of  the  seed  and  pump  beams,  Raman  gain, 
gain  saturation/pump  depletion,  RIVED,  and  nonlinear  dispersion. 
All  these  effects  give  rise  to  nonreciprocal  aberrations  in  the 
RA.  The  degree  to  which  a  phase  conjugate  reflector  can  compen¬ 
sate  for  reciprocal  aberrations,  in  the  presence  of  nonreciprocal 
aberrations,  has  been  evaluated  and  was  discussed  in  detail  in 
the  preceding  paragraphs.  The  results  of  this  study  can  be 
summarized  as  follows: 
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(1)  We  modeled  nonreciprocal  phase  errors  present  in  RAs  and 
found  very  few  situations  in  which  the  phase  conjugator 
was  not  able  to  compensate  for  simultaneously  induced 
reciprocal  phase  errors.  In  these  few  situations,  we 
foun-d  that  the  Strehl  ratio  of  the  conjugated  beam  would 
drop  by  a  factor  of  2  or  3  from  what  it  would  be  if 
there  were  no  reciprocal  aberrations  in  the  system.  In 
these  cases  we  are,  however,  modeling  reciprocal 
aberrations  strong  enough  to  reduce  the  Strehl  ratio  by 
at  least  2  to  3  orders  of  magnitude  in  a  system  without 
a  conjugator. 

(2)  When  phase  aberrations  of  the  seed  beam  are  strong 
enough  to  distort  its  intensity  profile  prior  to  con¬ 
jugation,  preservation  of  amplitude  information  during 
conjugation  is  required.  High  gain  in  the  RA  can  lead 
to  a  strongly  depleted  pump  before  significant  diffrac¬ 
tion  has  occurred,  and  has  the  same  effect  on  beam 
quality  as  a  conjugator  that  does  not  preserve  amplitude 
information . 

(3)  Pump-dependent  refractive  indexes  impair  the  ability  of 
the  conjugator  to  compensate  for  seed  aberrations. 
Dispersion  from  operation  off  line-center  of  the  Raman 
gain  does  not  produce  serious  effects.  However,  RIVED, 
resulting  from  the  transfer  of  significant  numbers  of 
molecules  to  their  excited  states,  can  be  important  for 
short  pulses.  This  requires  intense  pulses  shorter  than 
the  relaxation  time  of  the  excited  vibrational  level  and 
the  transit  time  of  the  system,  and  becomes  noticeable 
for  pulse  energy  densities  >~1  J/cm2 . 

(4)  For  high  gain  and  low  Fresnel  number  geometries,  strong 
aberrations  can  occur  from  nonuniform  gain.  Any 
impairment  of  the  ability  of  the  conjugation  process  to 
compensate  for  reciprocal  phase  errors  was  masked  by 
deterioration  of  the  beam  due  to  nonuniform  gain. 


SECTION  4 


EXPERIMENTAL  RESULTS 

4.1  GOALS  OF  THE  EXPERIMENTAL  PROGRAM 

The  goal  of  the  experimental  effort  on  this  program  was  to 
set  up  a  double-pass  RA/SBS  conjugator  system,  and  evaluate  its 
performance.  More  specifically,  the  objectives  were  to: 

(1)  demonstrate  that  SBS  phase  conjugation  can  be  used  to  compen¬ 
sate  for  reciprocal  phase  errors  in  the  double-pass  RA, 

(2)  demonstrate  that  nonreciprocal  phase  errors  exist  in  these 
systems,  (3)  evaluate  the  effect  of  nonreciprocal  phase  errors  on 
system  performance,  and  (4)  compare  the  experimental  results  with 
the  simulation  code  described  in  the  previous  section. 

In  previous  programs  at  the  Hughes  Research  Laboratories 
(HRL) 1 3 ’ 2 5 • 2 6 ’ 2 7  and  in  the  published  literature,28  SBS  has  been 
shown  to  produce  high  fidelity  conjugate  reflections.  Unfortun¬ 
ately,  quantitative  results  and  evaluation  of  the  range  of  valid¬ 
ity  are  not  easily  obtained.  A  major  effort  during  the  experi¬ 
mental  program  was  devoted  to  quantifying  the  beam  fidelity  so 
that  the  effect  of  nonreciprocal  phase  errors  could  be  demon¬ 
strated.  After  arriving  at  methods  to  quantify  the  beam 
fidelity,  we  proceeded  to  demonstrate  the  dependence  of  the 
fidelity  on  the  RA  characteristics. 

4.2  DATA  REDUCTION  AND  INTERPRETATION  OF  RESULTS 

Our  experimental  program  consisted  of  measuring  the  seed  beam 
fidelity  as  a  function  of  experimentally  adjustable  parameters 
for  both  multimode  and  single-mode  laser  operation.  The  adjust¬ 
able  parameters  include  the  pump  energy,  the  seed  energy,  and  the 
amount  of  pump  and  seed  aberration.  The  dependence  of  the 
fidelity  on  these  parameters  gives  insight  into  the  physical 
processes  that  limit  the  fidelity  in  the  double-pass  RA/SBS 
system.  The  beam  fidelity  was  determined  from  near-  and  far- 
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field  energy  distributions  of  the  pump  and  seed  beams.  These 
distributions  were  measured  using  the  laser  intensity  profiling 
system  (LIPS)  discussed  below. 

The  LIPS  used  during  these  experiments  consisted  of  a  photo 
diode  array  camera,  video  digitizer,  video  frame  buffers,  and  a 
computer  for  synchronization  and  data  analysis.  Figure  18  is  a 
schematic  of  the  LIPS  system.  Optical  beams  are  recorded  using 
an  EG&G  Reticon  camera  (m/n  MC9256) ,  which  has  a  256x256  photo¬ 
diode  array,  controlled  by  an  EG&G  formatter  (m/n  MB900)  to  pro 
duce  an  RS-170  compatible  video  signal.  The  video  signal  is 
digitized  by  an  Imaging  Technology  analog  processor  (m/n  AP512) 
and  stored  in  one  of  three  Imaging  Technology  frame  buffers 
(m/n  FB512) .  Each  frame  buffer  is  capable  of  holding  a  single 
512x512  video  image,  or  up  to  four  256x256  video  images.  A  com 
panion  Imaging  Technology  arithmetic  logic  unit  (m/n  ALU512) 
assists  the  analog  processor  for  high  speed  video  manipulations 
The  video  digitizer  and  optical  experiment  are  controlled  by  a 
Digital  Equipment  Corporation  PDP  micro-11/73  computer  with  a 
30  Mbyte  hard  disk,  a  dual  5-1/4-in.  floppy  drive,  a  streamer 
tape  unit,  and  a  LA50  printer  and  VT240  graphics  terminal  for 
input  and  output.  Synchronization  for  the  camera  and  the 
computer 
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Figure  18.  Schematic  of  the  laser  intensity  profiling 
system  used  to  evaluate  near-  and  far-field 
intensity  patterns. 


system  is  provided  by  the  video  digitizer.  The  computer  system 
synchronizes  the  pulsed  lasers  to  the  video  system.  Extensive 
software  has  been  written  on  this  system  for  controlling  optical 
experiments  and  data  reduction. 

The  LIPS  digitizes  the  energy  density,  for  a  single  shot  of 
the  laser,  on  a  256x256  grid.  Several  measures  of  beam  quality 
are  possible,  including  near-  and  far-field  intensity  profiles, 
beam  sizes,  Strehl  ratios,  and  times-di f f raction-1 imited  (XDL) . 
Near-  and  far-field  intensity  profiles  are  the  basic  data  from 
which  we  calculate  the  beam  fidelity.  They  are  output  either  as 
contour  plots  of  the  intensity,  or  as  profiles  taken  in  the  hori¬ 
zontal  and/or  vertical  directions.  Typical  far-field  beam  con¬ 
tour  and  profile  plots  are  shown  in  Figures  19  and  20. 

These  plots  give  good  qualitative  indication  of  the  beam 
shape,  but  not  a  quantitative  measure.  We  have  used  a  "profile 
index"  as  a  quantitative  measure  of  the  beam  shape.  This  index, 

0  ,  is  defined  by 

n  =  «ixi3>  -  <i2>2)/<i2>2  ,  (17) 

where  I  is  the  energy  density  and  the  brackets  "<>"  indicate  an 
integration  over  the  transverse  coordinates.  This  profile  index 
is  zero  for  a  flat-top  beam,  is  one-third  for  a  Gaussian  beam, 
and  is  0.490  for  the  Airy  disk  pattern.  Generally,  the  more 
spiked  the  beam  is,  the  higher  the  profile  index.  Figure  21 
shows  a  bullet  plot  of  the  profile  indexes  measured  for  the  beams 
analyzed  under  this  program.  The  ticks  labeled  "G"  and  "A"  indi¬ 
cate  the  profile  index  of  a  Gaussian  beam  and  an  Airy  disk  pat¬ 
tern,  respectively.  The  solid  and  open  triangles  represent  the 
near-  and  far-field  multimode  laser  measurements,  respectively. 
The  solid  and  open  diamonds  represent  the  near-  and  far-field 
single-mode  laser  experiments.  As  can  be  seen  from  the  figure, 
the  profile  index  of  the  near-field  exposures  lie  between  the 
Gaussian  and  the  flat-top  results.  The  multimode  far-field 
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F igure  21 . 


Bullet  plot  of  the  profile  index  for 

for  the  beams  analyzed  during  this  program. 

A  flat-top  beam  has  a  profile  index  of  zero, 
a  Gaussian  beam  has  an  index  of  one  third 
(indicated  by  ’G’),  and  an  Airy  disk  an  index 
of  0.490  (indicated  by  ’A’) . 


results  lie  generally  between  the  Gaussian  and  Airy  index,  and 
the  single-mode  far-field  index  is  generally  higher.  This 
implies  that  our  beams  are  not  strictly  Gaussian  or  flat-top  but 
actually  lie  somewhere  between. 

The  fidelity  of  the  laser  beam  was  measured  via  the  Strehl 
ratio  and  the  XDL .  Both  of  these  measures  compare  the  observed 
far-field  to  the  far-field  of  an  equivalent  unaberrated  beam. 

For  flat-top  and  Gaussian  near-fields,  simple  solutions  exist  for 
the  far-field  beam  shapes.  However,  reference  to  Figure  21  indi¬ 
cates  that  our  beams  are  neither  Gaussian  nor  flat-top,  and  the 
assignment  of  an  equivalent  far-field  beam  shape  is  difficult. 

In  addition,  this  assumption  dramatically  affects  the  calculated 
far-field  of  the  beam;  that  is,  the  far-field  radius  of  the 
Gaussian  beam  is  approximately  half  as  large,  while  its  on-axis 
intensity  is  four  times  larger.  To  make  our  results  comparable 
to  the  work  of  others,  we  have  evaluated  the  beam  fidelity  both 
under  assumed  near-field  shapes  for  the  multimode  laser  results, 
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and  under  the  Gaussian  near-field  shape  for  the  single-mode 
resu 1 tr  . 

A  beam  with  a  flat-top  near-field  diffracts  into  an  Airy  disk 
pattern  in  the  far-field.  When  examined,  by  observing  in  the 
focal  plane  of  a  lens  with  focal  length  f,  the  far-field  inten¬ 
sity  pattern  is  given  by 

Ijf(p)  =  4l?f  J2 [(kp?f/f)p]/ [(kp?f/f)p] 2  ,  (18) 

where  subscripts  indicate  the  flat-top  beam,  superscripts  the 
far-field  (ff)  and  near-field  (nf),  I£f  intensity  in  the  near¬ 
field,  Jj  a  first  order  Bessel  function,  k  =  27T/X,  and  the 

near-field  radius  of  the  beam.  We  model  the  near-field  radius  of 
our  "flat-top"  beam  by  setting  p£f  equal  to  the  radius  of  an 
aperture  that  contains  90%  of  the  energy  in  our  near-field 
exposure.  The  far-field  "flat-top,"  pff,  radius  is  taken  as  the 
aperture  radius  that  contains  83.79%  of  the  energy  of  the  far- 
field  exposure,  the  same  fraction  of  the  energy  as  contained 
within  the  first  dark  ring  of  the  Airy  disk  pattern.  In  both 
cases  the  aperture  is  positioned  to  allow  maximum  throughput  from 
a  smaller  radius  giving  ~50%  transmission.  To  obtain  the  XDL , 
the  measured  far-field  radius  (p£f)  is  compared  t,o  the  radius  of 
the  first  zero  in  the  Airy  disk  pattern  calculated  from  the  near¬ 
field  radius  (£>ff).  That  is, 

XDLf  -  pJVCii.i  fAP?f)  ,  (19) 

where  j  is  the  first  positive  zero  of  the  first  order  Bessel 
function  (j l  j  ~3. 83171).  The  Strehl  ratio  is  calculated  by 
comparing  the  highest  intensity  pixel  in  the  exposure  with  the 
peak  intensity  calculated  from  the  Airy  disk  pattern,  both  nor¬ 
malized  to  the  total  power  in  the  beam.  It  follows  that  the 
Strehl  ratio  8f  is  given  by 


0f  -  (a/E0)/[7r(p?f)2/2X2f2] 
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where  a  is  the  exposure  level  in  the  hottest  Dixel,  and  E0  is  the 
total  exposure  of  the  beam. 

A  beam  with  a  Gaussian  near-field  diffracts  into  another 
Gaussian  in  the  far-field.  When  examined,  by  observing  in  the 
focal  plane  of  a  lens  with  focal  length  f,  the  far-field  inten¬ 
sity  pattern  is  given  by 


I V  (p)  --  I"f(0)  >(p"f)  Z/Xf  ]  2  exp  (  -  2  [9/3/  (Xf/?r/7"T)]  2  )  ,  (21) 


n  f 


,  n  f  \  2 


>-1  2 


n  f  ' 


where  the  subscripts  indicate  the  Gaussian  beam.  We  model  the 
radii  of  the  ’’Gaussian"  beams  (p£f  and  p **)  by  setting  them  equal 
to  the  radii  of  an  aperture  that  contains  86.47%  of  the  energy  in 
the  exposures,  the  same  fraction  of  the  energy  as  contained 
within  the  1/e2  radius  of  a  perfect  Gaussian  beam.  As  above,  the 
aperture  is  positioned  to  allow  maximum  throughput  from  a  smaller 
radius  giving  ~50%  transmission.  To  obtain  the  XDL ,  the  measured 


far-field  radius  ( p  9J  is  compared  to  the  1/e2  radius  of  the  far- 


field  Gaussian  beam  calculated  from  the  near-field  radius  (p £  ). 


That  ii 


XDLq  =  plf/(Xf/*pr) 


n  f  • 


(22) 


The  Strehl  ratio  is  calculated  by  comparing  the  highest 
intensity  pixel  in  the  exposure  with  the  peak  intensity  cal¬ 
culated  from  the  far-field  Gaussian  beam  pattern,  both  normalized 
to  the  total  power  in  the  beam.  It  follows  that  the  Strehl  ratio 


8q  is  given  by 


9  a  =  (a/E0)  /  [27r(p"f/\f)  2  ] 


(23) 


Our  experimental  program  consisted  of  evaluating  the  profile 
indexes,  the  near-  and  far-field  radii,  the  Strehl  ratios,  and 
the  XDLs  as  a  function  of  the  experimentally  adjustable  para¬ 
meters.  Comparison  with  theoretical  and  simulation  results 
allows  identification  of  the  nonreciprocal  aberrations. 
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4.3  MULTIMODE  LASER  RESULTS 


Multimode  laser  investigations  were  carried  out  using  the 
experimental  apparatus  depicted  in  Figure  22.  A  YG581C  YAG  laser 
manufactured  by  Quantel  International  was  used  as  the  optical 
source.  This  laser  consists  of  an  actively  Q-switched  NdtYAG 
oscillator/amplifier,  designed  for  high  energy  applications. 

This  results  in  an  energetic,  but  highly  multimode  (transverse 
and  longitudinal)  laser  output.  The  ~700  mJ  IR  (X  =  1.06  /xm) 
output  of  the  oscillator/amplifier  is  frequency  doubled,  result¬ 
ing  in  ~300  mJ  of  vertically  polarized  optical  energy  in  the 
visible  spectrum  (X  =  532  nm) .  This  output  is  split  by  a  30% 
reflecting  beam  splitter;  the  70%  transmitted  is  used  to  pump 
backward  stimulated  Raman  scattering  (BSRS)  for  generating  the 
seed  beam  (see  below) ,  and  the  30%  reflected  by  the  beam  splitter 
is  used  to  pump  the  RA . 

The  polarization  of  the  optical  pulse  transmitted  by  the  beam 
splitter  is  first  rotated  to  the  horizontal  plane  by  a  half-wave 
plate  (X/2),  and  then  focused  into  a  high  pressure  H2  cell  to 
generate  the  Raman  seed  beam  by  BSRS.  This  red  (X  =  683  nm)  seed 
beam  is  separated  from  the  green  pump  beam  using  an  optical 
prism.  The  seed  beam  then  passes  through  a  3:5  telescope  reducer 
and  an  optional  aberration  plate  before  entering  the  RA .  During 
its  first  pass  through  the  RA ,  the  seed  beam  senses  aberrations 
due  to  optics  and  inhomogeni ties  in  the  amplifier.  Because  the 
seed  beam  is  polarized  orthogonal  to  the  pump  beam  on  the  first 
pass,  the  Raman  gain  is  very  small;  in  addition,  the  Raman  gain 
could  be  zero  because  the  pump  beam  has  not  entered  the  ampli¬ 
fier.  After  leaving  the  amplifier,  the  seed  passes  through  a 
dichroic  mirror  and  is  converted  to  circular  polarization  by  a 
quarter-wave  plate  (X/4) .  The  resulting  beam  is  focused  into  a 
SBS  phase  conjugation  cell  and  reemerges  with  circular  polariza¬ 
tion  of  the  opposite  handedness.  Thus,  the  polarization  of  this 
conjugated  beam  is  converted  to  vertical  linear  polarization  by 
the  quarter-wave  plate  prior  to  being  combined  with  the  pump  at 


YAG  LASER 
=  532  nm 


3 


the  dichroic  mirror.  After  being  combined,  the  conjugated  seed 
and  the  coincident  pump  beam  enter  the  RA ,  transferring  energy 
from  the  pump  beam  to  the  seed  beam.  The  amplified  conjugated 
seed  beam  repasses  the  optional  aberration  plate  and  is  split  off 
in  a  polarizer  by  virtue  of  its  rotated  polarization.  At  this 
point,  the  beams  (seed  and/or  pump)  are  examined  to  determine 
their  near-  and  far-field  characteristics.  The  near-field  is 
measured  at  the  polarization  beam  splitter,  while  the  far-field 
is  measured  at  the  focal  plane  of  a  100  cm  lens. 

The  pump  beam  for  the  RA  is  split  off  by  the  30%  beam  split¬ 
ter  and  undergoes  an  optical  delay  by  transitting  the  table 
several  times.  The  optical  delay  of  the  pump  beam  is  adjusted  so 
the  pump  will  be  coincident  with  the  seed  beam  on  its  second  pass 
through  the  RA .  Two  pumping  configurations  are  used.  In  the 
first  conf iguration ,  a  long  focal  length  lens  is  placed  in  the 
optical  delay  path.  This  lens  images  the  laser  output  into  the 
RA,  producing  a  more  uniform  pump  intensity  than-  when  the  pump 
beam  is  allowed  to  simply  diffract  during  the  optical  delay.  In 
the  second  configuration,  the  pump  beam  is  allowed  to  diffract 
during  the  optical  delay. 

The  Quantel  laser  produces  a  highly  multimode  (transverse  and 
longitudinal)  optical  pulse  of  ~10  ns  duration.  During  these 
experiments,  the  laser  was  operated  in  two  different  modes: 
"smooth  pulse,"  and  "narrow  linewidth."  To  operate  this  laser  in 
the  smooth  pulse  mode,  an  intracavity  reflector  is  added  to  the 
oscillator.  The  effect  of  this  reflector  is  to  fill  in  the  mode 
spectrum  of  the  laser,  driving  the  intermode  beat  frequency 
beyond  the  bandwidth  of  common  laboratory  detectors.  The  result 
is  a  pulse  that  appears  temporally  "smooth."  This  reflector 
increases  the  number  of  oscillating  modes,  but  the  laser  band¬ 
width  remains  at  the  norma]  operation  level  of  ~30  GHz.  In  the 
narrow  linewidth  mode  of  operation,  an  etalon  is  added  to  the 
oscillator  to  reduce  the  mode  spectrum.  With  this  etalon,  the 
laser  bandwidth  is  reduced  to  ~7  GHz.  As  discussed  below,  we 
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found,  as  expected,  that  in  the  smooth  pulse  mode  the  efficiency 
of  the  stimulated  scattering  processes  was  reduced  because  of  the 
wide  bandwidth.  In  fact,  we  were  unable  to  obtain  sufficient 
reflection  from  the  SBS  conjugator  to  demonstrate  the  operation 
of  the  RA/SBS  system.  In  the  narrow  linewidth  mode  of  operation, 
the  laser  linewidth  was  sufficiently  narrow  to  obtain  good  SBS, 
and  operation  of  the  RA/SBS  system  was  possible. 

The  seed  beam  for  our  RA/SBS  system  was  created  by  BSRS . 

This  process  has  been  used  successfully  to  efficiently  generate 
high  quality  "conjugate"  seed  beams  using  narrow  bandwidth 
lasers.29  With  wide  bandwidth  pump  lasers,  the  BSRS  process  is 
less  efficient  and  Raman  seed  beams  are  typically  generated  using 
forward  stimulated  Raman  scattering  (FSRS)  in  single-  or  multi¬ 
pass  cells,  or  using  resonator  configurations.30  In  this  pro¬ 
gram  the  BSRS  and  the  FSRS  processes  were  investigated,  and  the 
results  are  shown  in  Figures  23  and  24.  In  Figure  23,  the  stimu¬ 
lated  scattering  output  is  shown  as  a  function  of  the  pumping  F- 
number.  At  low  F-numbers,  the  BSRS  output  dominates,  maximizing 
at  ~F/15  with  a  quantum  efficiency  of  ~17%.  and  beam  quality 
matching  the  pump.  At  high  F-numbers,  the  FSRS  output  dominates, 
achieving  a  quantum  efficiency  of  >60%  at  ~F/100.  The  beam 
quality  is,  however,  very  poor.  In  Figure  24,  the  stimulated 
Raman  scattering  (SRS)  output  is  shown  as  a  function  of  the  cell 
pressure  in  low  and  high  F-number  geometries.  In  the  low  F-num- 
ber  (F/13.6)  geometry,  the  BSRS  output  dominates  for  low  pressure 
and  is  in  strong  competition  with  the  FSRS  at  high  pressures.  In 
the  high  F-number  (F/104)  geometry,  the  FSRS  output  continues  to 
increase  with  pressure . 

Simultaneous  experiments,  discussed  below,  demonstrated  that, 
while  the  BSRS  process  worked  with  the  broad  bandwidth  laser,  the 
SBS  process  would  not.  Because  of  this,  the  line  narrowing 
etalon  was  inserted  into  the  laser,  narrowing  the  bandwidth  to 
~7  GHz.  With  the  new  narrower  linewidth,  the  BSRS  process  became 
much  more  efficient.  These  results  are  shown  in  Figures  25  and 
26.  In  Figure  25  the  BSRS  output  energy  is  shown  as  a  function 
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Energy  contained  in  an  SRS  generated  seed 
pulse.  The  hydrogen  filled  Raman  cell  was 
pumped  at  a  532  nm  wavelength  with  a  200  mJ 
wide  bandwidth  (30  GHz)  laser  pulse.  The 
circles  and  squares  show  the  BSRS  and  FSRS 
output,  respectively,  in  an  F/13.6  geometry. 
The  triangles  show  the  FSRS  output  in  an  F/104 
geometry . 
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Energy  contained  in  a  BSRS  generated  seed 
pulse.  The  hydrogen  Raman  cell  (1580  psi)  was 
pumped  at  a  532  nm  wavelength  with  a  narrow 
bandwidth  (7  GHz)  laser  pulse.  The  circles, 
squares,  and  triangles  show  the  BSRS  output  in 
an  an  F/12.5,  and  an  F/25  geometry, 

respectively. 
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Figure  26. 


Energy  contained  in  the  BSRS  generated  seed 
pulse.  The  hydrogen  filled  Raman  cell  was 
pumped  at  a  532  nm  wavelength  with  a  200  mJ 
narrow  bandwidth  (7  GHz)  laser  pulse  in  an 
F/12.5  geometry. 


of  the  pump  energy.  The  figure  shows  that  the  BSRS  process  is 
much  more  efficient  with  the  narrow  bandwidth  laser,  attaining  a 
quantum  efficiency  of  ~50%  in  an  F/12.5  configuration.  This  is 
almost  three  times  the  value  attained  with  the  wide  bandwidth 
laser.  Figure  26  shows  the  BSRS  output  as  a  function  of  cell 
pressure  for  a  F/12.5  pumping  geometry.  The  output  energy 
increases  with  cell  pressure,  saturating  at  ~1200  psi .  For  the 
purposes  of  demonstrating  nonreciprocal  aberrations  in  the  RA,  we 
operated  the  BSRS  seed  generator  at  a  pressure  of  ~2000  psi  with 
a  narrow  bandwidth  pump  of  ~200  mJ .  This  gave  a  seed  pulse  of 
~90  mJ  at  the  input  of  the  RA ,  and  was  large  enough  to  obtain 
sufficient  reflectivity  from  the  SBS  con j ugator . 

SBS  phase  conjugation  of  the  multimode  laser  beams  was  per¬ 
formed  in  a  high  pressure  gas  cell.  Both  methane  (CH4)  and  sul¬ 
fur  hexafluoride  (SFg)  were  evaluated  as  conjugator  materials. 
Neither  gave  suitable  reflectivities  with  the  wide  bandwidth 
laser.  In  fact,  the  SBS  reflection  was  sporadic  and  of  poor 
quality.  This  held  true  under  all  pressures,  focusing  condi¬ 
tions,  and  pump  energies  at  both  532  and  683  nm  pump  wavelengths. 

With  the  narrow  bandwidth  laser,  the  situation  was  consider¬ 
ably  improved.  An  SBS  phase  conjugate  signal  could  be  observed 
from  the  pump  beam  and  the  BSRS  seed  beam  in  both  CH4  and  SFg. 

The  observed  reflectivity  was  in  the  20%  to  30%  range  for  pump 
energies  ranging  from  25  to  200  m J .  At  any  given  pressure,  the 
SBS  gain  of  SFe  is  substantially  higher  than  CH4 .  However,  SF6 
liquefies  at  a  pressure  of  ~300  psi  (room  temperature),  and  at 
this  pressure,  the  SBS  gain  in  SFg  is  comparable  to  the  gain  in 
CH4  at  high  pressure  (~2000  psi).  In  initial  experiments  with 
SF 6  (300  psi),  we  observed  BSRS  in  addition  to  the  SBS  signal; 
undoubtly  the  presence  of  this  competing  nonlinear  process 
reduces  the  fidelity  of  the  phase  conjugate  signal.  No  BSRS  was 
observed  with  CH4 .  For  these  reasons,  we  selected  CH4  as  the 
best  material  for  the  SBS  phase  conjugator.  The  SBS  cell  was 
37  cm  long,  contained  CH4  at  a  pressure  of  2000  psi,  and  was 


operated  in  an  F/10  geometry.  The  SBS  reflectivity  of  our  CH4 
cell  was  estimated  by  first  measuring  the  system  losses  and  then 
taking  them  out  of  the  reflectivities  measured  for  the  data  dis¬ 
cussed  below.  These  estimates  are  plotted  in  Figure  27  as  a 
function  of  input  to  the  SBS  cell.  As  can  be  seen,  the  SBS 
reflectivity  is  typically  in  the  range  of  20  to  30 %  for  input 
energies  of  >20  m J .  At  higher  input  energies,  the  SBS 
reflectivity  did  not  increase  substantially. 

The  RA  used  in  these  experiments  was  a  107-cm-long  high  pres¬ 
sure  cell  containing  -2000  psi  of  hydrogen.  At  this  pressure, 
the  Raman  lineshift  A u r  =  4154  80  cm'1  gives  a  peak  Raman  gain  at 
683  nm  (pump  at  532  nm) .  The  Raman  linewidth  Sur  is  approxi¬ 
mately  7.1  GHz  at  this  pressure.31  The  SBS  frequency  s.-^ft  for 
our  CH4  cell  is  only  1.3  GHz  and  the  effect  of  being  off  line- 
center  of  the  Raman  gain  is  negligible.  The  peak  Raman  gain 
gG  =  3.0xl0'9  cm/W  gives  a  small-signal  gain  coefficient 
g0 IL  -2.4  f  or  a  40  mJ  pump  pulse. 

As  discussed  above,  our  experimental  observation  of  nonrecip¬ 
rocal  aberrations  consisted  of  evaluating  the  amplified  seed 
fidelity  as  a  function  of  experimental  parameters.  We  begin  this 
discussion  by  characterizing  the  pump  beam,  the  seed  beam,  and 
the  conjugation  process.  This  is  followed  by  a  discussion  of  the 
effect  of  Raman  gain  and  saturation  of  Raman  gain  on  the  conju¬ 
gate  seed  beam  fidelity. 

As  discussed  above,  the  fidelity  of  our  beams  was  characteri¬ 
zed  by  evaluating  the  size  of  the  beam  and  its  peak  intensity. 

The  results  of  these  evaluations  on  the  pump  beam  are  shown  in 
Figures  28,  29,  and  30.  In  Figure  28,  the  near-  and  far-field 
pump  beam  radii  are  shown  as  a  function  of  the  pump  energy. 

These  beam  profiles  were  evaluated  after  the  pump  beam  traversed 
the  RA  in  the  absence  of  the  seed  wave.  As  such,  variations  in 
the  pump  beam  can  be  attributed  to  nonlinear  interaction  of  the 
pump  with  itself  (e  g.  ,  se 1 f - focus i ng)  .  The  radii  are  given  by 
the  radius  of  an  aperture  that,  when  centered  on  the  beam,  trans¬ 
mits  a  fixed  fraction  of  the  beam.  The  0.9000  (triangles)  and 
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Figure  27. 


Estimated  reflectivity  of  CH4  SBS  cell 

at  a  pressure  of  2000  psi  for  narrow  bandwidth 

(~7  GHz)  radiation  at  a  wavelength  of  683  nm. 
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Figure  28 


Pump  beam  radius  as  a  function  of  the  energy 
of  the  pump  pulse.  The  radii  are  given  by  the 
radius  of  an  aperture  that,  when  centered  on  the 
beam,  transmits  a  fixed  fraction  of  the  beam. 

The  triangles  and  diamonds  give  the  radii 
corresponding  to  0.9000  and  0.8647  of  the  near- 
field  beam,  while  the  circles  and  squares  give 
radii  corresponding  to  0.8379  and  0.8647  of  the 
far-field  beam,  respectively.  The  solid  symbols 
are  for  the  imaged  pumping  geometry,  while  the 
open  symbols  are  for  the  diffracted  pumping 
geometry . 


the  0.8647  (diamonds)  radii  of  the  near-field  beam  are  used  to 
estimate  the  radius  of  an  equivalent  flat-top  and  Gaussian  beam, 
respectively.  The  solid  symbols  in  the  figure  indicate  results 
obtained  in  the  "imaged"  pumping  geometry,  and  the  open  symbols, 
the  "diffracted"  pumping  geometry.  The  value  for  a  pump  energy 
of  "0"  was  obtained  by  removing  the  Raman  cell .  The  data  show 
very  little  change  in  the  near-field  beam  radii.  The  small 
increase  in  the  beam  radii  observed  for  high  pump  energy  indi 
cates  that  a  nonlinear  mechanism,  such  as  se 1 f - de focus i ng ,  may  be 
present  in  the  amplifier,  although  its  effect  is  small.  The 
near-field  radius  of  the  pump  beam  was  estimated  by  a  linear 
least  squares  extrapolation  to  a  pump  energy  of  0.0  mJ .  These 
values  were  used  to  estimate  the  far-field  radii  and  peak  inten¬ 
sities  as  listed  in  the  table  below. 

TABLE  1.  Characteristics  of  Multimode  Pump  Beam. 


PUMP  BEAM 

IMAGED 

GAUSSIAN  |  FLAT  TOP 

DIFFRACTED 
GAUSSIAN  |  FLAT  TOP 

NEAR  FIELD 
RADIUS  (cm) 

r 

0.419 

0.440 

0  504 

0  530 

FAR  FIELD 
RADIUS  (cm) 

0  00404 

0.00738 

0  003.  :6 

i_  ..... 

0.00612 

Iff (0)/Eo 
( cm  1 ) 

38919 

10739 

— 

56334 

15587 

The  XDL  of  the  pump  beam  is  plotted  as  a  function  of  the  pump 
energy  in  Figure  29.  The  circles  show  the  flat  top  results  and 
the  squares  the  Gaussian  results.  As  above,  the  open  symbols 
represent  the  diffracted  pumping  configuration,  and  the  solid 
symbols  represent  the  imaged  pumping  configuration  As  expected, 
the  XDL  of  the  pump  beams  increases  slightly  with  increasing  pump 
energy.  In  Figure  30,  the  Strehl  ratio  of  the  pump  beam  is  plot 
ted  as  a  function  of  the  pump  energy  The  Strehl  ratio  is  seen 
to  initially  increase  and  then  to  decrease  This  amounts  to  an 


~2X  change  in  the  Strehl  ratio.  Reference  to  Figure  29  indicates 
that  the  far-field  radius  (XDL)  changes  by  <20%,  accounting  for 
an  ~40%  change  in  the  Strehl  [Strehl  a  1/(XDL)2] .  Because  the 
change  in  Strehl  is  larger,  the  far-field  intensity  profile  must 
also  be  changing.  That  is,  at  high  pump  energies  the  far-field 
energy  distribution  flattens  out. 

The  quality  of  our  seed  beam  was  characterized  in  an  equi¬ 
valent  manner  to  the  pump  beam  characterization  discussed  above. 
This  characterization  was  done  without  the  RA ,  the  pump  beam,  or 
the  SBS  phase  conjugator;  the  results  are  given  below: 


TABLE  2.  Characteristics  of  Multimode 
Seed  Beam. 


SEED 

BEAM 

GAUSSIAN 

CALCULATION 

FLAT-TOP 

CALCULATION 

NEAR  FIELD 
RADIUS  (cm) 

0.406 

0 . 428 

FAR  FIELD 
RADIUS  (cm) 

0.00536 

0.00974 

Iff  (0)/Eo 

(cm'2) 

22160 

6163 

<XDL> 

23.73 

12.46 

< STREHL > 

0.00206 

0.00740 

The  results  labeled  <XDL>  and  <STREHL>  are  the  XDL  and  Strehl 
ratio  calculated  for  12  intensity  profiles  taken  without  the  RA . 
These  results  are  used  in  all  subsequent  evaluations  of  the  seed 
beam  fidelity.  In  this  manner,  the  seed  fidelity  is  calculated 
relative  to  that  observed  without  a  RA  or  SBS  phase  conjugator. 

In  a  similar  manner,  the  fidelity  of  the  phase  conjugate  beam  is 
evaluated.  There  was  no  significant  difference  between  these 
results  and  those  given  above.  To  demonstrate  that  the  SBS  phase 
conjugator  can  compensate  for  reciprocal  phase  errors  in  the  RA , 
an  aberration  plate  was  inserted  between  the  dichroic  mirror  and 


the  quarter  wave  plate.  This  location  was  chosen  because  the 
plate  induced  such  large  aberrations  that  the  seed  beam  would  be 
clipped  by  the  optics  if  the  aberrator  had  been  placed  further 
away  from  the  SBS  cell  focusing  lens.  With  this  aberrator,  the 
Strehl  ratio  of  the  seed  beam  reflected  from  an  ordinary  mirror 
decreased : 

TABLE  3.  Characteristics  of  Aberrated/ 

Unconjugated  Multimode  Seed 
Beam . 


ABERRATED 
SEED  BEAM 

GAUSSIAN 

CALCULATION 

FLAT-TOP 

CALCULATION 

NEAR  FIELD 
RADIUS  (cm) 

0.557 

0.581 

FAR-FIELD 
RADIUS  (cm) 

0.387 

0.367 

<XDL> 

00 

rH 

IS 

37.75 

<STREHL> 

0 .00025 

0.00091 

When  reflected  from  the  SBS  phase  conjugator,  the  fidelity  of  th 
seed  beam  was  much  improved  (see  Table  4) . 

Comparison  of  the  data  in  these  three  tables  indicates  that 
while  the  phase  conjugate  fidelity  was  not  perfect,  the  phase 
conjugator  was  able  to  compensate  for  the  most  of  the  reciprocal 
phase  error  introduced  by  the  aberration  plate.  In  the  discus¬ 
sions  that  follow,  we  observe  reductions  in  the  Strehl  of  <50%; 
we  believe  the  SBS  conjugator  is  capable  of  providing  a  good 
conjugate  of  beams  with  aberrations  of  this  size.  In  addition, 
we  investigate  the  change  in  beam  fidelity  as  a  function  of  gain 
in  the  RA ;  as  a  result,  the  observed  changes  in  Strehl  must 
result  from  nonreciprocal  processes  in  the  amplifier. 

Our  analysis  of  nonreciprocal  aberrations  consisted  of  first 
evaluating  the  effect  of  Raman  gain,  and  second,  the  effect  of 
Raman  gain  saturation.  To  investigate  the  effects  of  Raman  gain 
*•  replaced  the  RA  cell  and  adjusted  the  pump  to  be  coincident 
•hr  conjugated  seed  beam.  In  these  experiments,  a  seed 


TABLE  4.  Characteristics  of  Aberrated/ 

Conjugated  Multimode  Seed  Beam. 


CONJUGATED 
SEED  BEAM 

GAUSSIAN 

CALCULATION 

FLAT-TOP 

CALCULATION 

NEAR  FIELD 
RADIUS  (cm) 

0.452 

0.476 

FAR-FIELD 
RADIUS  (cm) 

0 . 149 

0.140 

<XDL> 

27.86 

14.40 

<STREHL> 

0.00157 

0.00566 

pulse  of  ~30  mJ  entered  the  RA  on  the  first  pass.  Taking  the 
reflection  losses  of  the  RA  cell  windows,  the  dichroic  beam¬ 
splitter,  the  focusing  lens,  and  the  SBS  cell  windows,  we  esti¬ 
mate  that  ~17  mJ  of  seed  energy  actually  entered  the  SBS  cell. 
From  Figure  27,  the  anticipated  SBS  reflectivity  is  ~15%,  and, 
with  additional  system  losses,  we  estimate  that  ~1.5  mJ  of  con¬ 
jugated  seed  energy  actually  entered  the  RA  on  the  second  pass. 
Figure  31  shows  the  measured  Raman  gain  (Amplified  Stokes 
energy/input  seed  energy)  as  a  function  of  pump  energy.  At  the 
highest  gains,  we  estimate  that  ~30%  of  the  pump  energy  has  been 
depleted  by  the  amplification  process.  It  follows  that  pump 
depletion  is  a  consideration  in  these  measurements.  This  is  also 
evident  from  the  linear  increase  of  gain  with  pump  energy.  At 
very  low  gains,  pump  depletion  becomes  less  important  and  the 
gain  is  exponential  with  pump  energy.  Reference  to  the  figure 
indicates  that  this  is  true  only  for  pump  energies  <30  mJ. 

Figures  32  and  33  show  the  results  of  amplified  seed  beam 
fidelity.  In  Figure  32,  the  near-  and  far-field  seed  beam  radii 
are  shown  as  a  function  of  the  pump  energy.  These  beam  profiles 
were  evaluated  at  the  RA  output,  and  variations  in  the  radii  can 
be  attributed  to  aberrations  produced  by  the  interaction  of  the 
pump  and  seed  beams  within  the  RA .  The  near-field  radii, 
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PUMP  ENERGY,  mJ 


Seed  beam  radius  as  a  function  of  the 
pump  energy.  The  symbols  correspond  to 
used  in  Figure  28. 
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PUMP  ENERGY,  mJ 


Strehl  ratio  of  the  ampl i f ied-con j ugated  seed 
beam  as  a  function  of  the  pump  energy.  The 
symbols  correspond  to  those  used  in  Figure  30 


indicated  by  the  triangles  and  the  diamonds,  show  an  initial 
decrease  followed  by  an  increase.  The  initial  decrease  is 
thought  to  be  a  result  of  nonuniform  gain  caused  by  the  Gaussian- 
like  profile  of  the  pump  beam.  For  pump  energies  in  excess  of 
30  m J ,  the  Raman  gain  begins  to  saturate,  depleting  the  high 
intensity  center  of  the  pump  beam.  As  a  result,  near-field  radii 
increase  in  a  manner  much  like  that  discussed  in  Section  2.3. 

The  far-field  radii,  shown  by  the  circles  and  squares,  increase 
monotonical ly  with  pump  energy. 

In  Figure  34,  the  XDL  of  the  amplified  conjugated  seed  beam 
is  plotted  as  a  function  of  the  pump  energy.  As  is  evident  from 
Figure  34,  the  XDL  of  the  beam  shows  a  significant  (~50%) 
increase  with  Raman  gain.  Figure  33  shows  the  Strehl  ratio  of 
the  amp  1 i f i ed-con j ugated  seed  beam.  As  expected,  the  Strehl 
ratio  decreases  with  increasing  pump  energy.  The  initial  strong 
decrease  in  Strehl  is  a  result  of  nonuniform  gain,  as  discussed 
above.  At  pump  energies  above  ~30  m J ,  the  pump  beam  begins  to 
saturate  and  the  effects  of  nonuniform  gain  decrease.  At  this 
point,  the  Strehl  ratio  remains  roughly  constant  with  increasing 
energy.  At  the  highest  gains,  the  Strehl  ratio  of  the  amplified- 
conj ugated  seed  beam  (for  imaged  pumping  condition)  is  roughly 
equivalent  to  that  of  the  pump  beam,  because  the  pump  and  seed 
have  similar  profiles  under  these  conditions.  When  the 
diffracted  pumping  condition  is  used,  the  Strehl  ratios  drop  by 
nearly  a  factor  of  2;  this  drop  results  from  the  different  pump 
and  seed  profiles. 

To  investigate  the  effects  of  Raman  gain  saturation,  we 
evaluated  the  conjugate  beam  fidelity  as  a  function  of  seed 
energy.  In  these  experiments,  a  constant  pump  pulse  energy  of 
~40  mJ  entered  the  RA .  The  seed  energy  entering  the  RA  on  the 
second  pass  was  estimated  by  measuring  the  conjugated  seed  beam 
output  without  the  pump  beam  and  compensating  for  the  system 
losses.  Figure  35  shows  the  measured  Raman  gain  as  a  function  of 
seed  energy.  For  the  smallest  input  energy,  we  estimate  that  <1% 
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PUMP  ENERGY,  mJ 


XDL  of  the  amplified  conjugated  seed  beam 
as  a  function  of  the  pump  pulse  energy.  The 
symbols  correspond  to  those  used  in  Figure  30. 


SEED  ENERGY,  mJ 


Figure  35.  Measured  gain  in  the  RA  as  a  function  of  seed 
energy  entering  the  amplifier  on  the  second 
pass.  The  pump  energy  is  ~40  mJ.  The  symbol 
correspond  to  those  used  in  Figure  31. 
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of  the  pump  energy  has  been  transferred  to  the  seed  beam,  and 
thus  pump  depletion  plays  no  role.  At  the  highest  seed  energy, 
~60%  of  the  pump  energy  is  transferred  to  the  conjugated  seed, 
with  pump  depletion  playing  an  important  role.  Pump  depletion 
effect  on  conjugate  seed  fidelity  is  shown  in  the  next  three 
f igures . 

In  Figure  36,  the  near-  and  far-field  seed  beam  radii  are 
shown  as  a  function  of  the  seed  energy.  These  beam  profiles  were 
evaluated  at  the  RA  output,  and  variations  in  the  radii  can  be 
attributed  to  aberrations  produced  by  the  interaction  of  the  pump 
and  seed  beams  within  the  RA .  The  near-field  radii,  indicated  by 
the  triangles  and  the  diamonds,  show  a  gradual  increase  with  seed 
energy,  as  a  result  of  the  amplified  seed  beam  taking  on  a  more 
f 1  at- top- 1 i ke  profile  because  of  pump  depletion  in  the  hot  center 
of  the  seed.  This  change  in  profile  is  also  indicated  by  a 
reduction  in  the  profile  index  as  the  seed  energy  is  increased. 
The  far-field  radii,  shown  by  the  circles  and  squares,  also 
increase  monotonical ly  with  pump  energy. 

The  XDL  of  the  ampl  i  f  i ed-con  j  ugated  seed  beam  is  plotted  as  a 
function  of  the  seed  energy  in  Figure  37.  As  is  evident  from 
Figure  37,  the  XDL  of  the  beam  shows  an  increase  as  the  Raman 
gain  saturates.  The  observed  change  in  XDL  is  smaller  than  that 
observed  in  Figure  34,  which  is  probably  the  result  of  a 
decreased  conjugation  fidelity  for  low’  seed  energies 
(reflectivity  <20%,  for  seed  energy  <2  mJ) .  Figure  38  shows  the 
Strehl  ratio  of  the  amplif ied-conj ugated  seed  beam.  As  expected, 
the  Strehl  ratio  decreases  with  increasing  pump  energy.  The 
reduction  in  Strehl  ratio  shown  in  Figure  38  is  similar  to  that 
in  Figure  34,  and  is  believed  to  be  due  to  similar  processes. 

The  initial  strong  reduction  is  due  to  nonuniform  gain,  and  the 
subsequent  slow  reduction  is  due  to  pump  beam  saturation. 

In  summary,  we  have  observed  a  reduction  in  the  Strehl  ratio 
of  our  multimode  seed  beam,  as  a  result  of  amplification  by  a 
multimode  pump  beam  in  a  RA .  At  low  extraction,  when  pump 
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XDL  of  the  amplif ied-conj ugated  seed  beam 
as  a  function  of  the  seed  pulse  energy.  The 
symbols  correspond  to  those  used  in  Figure  29 
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Strehl  ratio  of  the  ampl i f ied-con j ugated  seed 
beam  as  a  function  of  the  seed  energy.  The 
symbols  correspond  to  those  used  in  Figure  30 


depletion  is  not  important,  large  decreases  in  the  beam  fidelity 
are  observed  as  a  result  of  nonuniform  gain,  in  turn  the  result 
of  a  pump  beam  that  is  not  strictly  flat-top.  As  the  pump 
becomes  depleted,  the  seed  beam  begins  to  take  on  the  intensity 
profile  of  the  pump  beam.  Further  reductions  in  amplified  seed 
beam  fidelity  are  not  significant. 

4.4  SINGLE-MODE  LASER  EXPERIMENTS 

After  conducting  multimode  experiments,  we  proceeded  to 
evaluate  the  RA/SBS  system  for  a  di f f ract ion- 1 imi ted ,  single¬ 
longitudinal  mode  (SLM)  pump  laser  with  a  532  nm  output  energy  of 
~250  m J .  For  these  experiments,  we  were  able  to  use  a  laser 
previously  developed  for  another  program.32  The  laser  consists 
of  a  commercial  Nd:YAG  laser  modified  for  operation  in  a  phase 
conjugate  master  oscillator  power  amplifier  (MOPA)  configuration 
as  shown  in  Figure  39.  The  oscillator  and  amplifier  heads  are 
from  a  QUANTEL  YG481  laser,  and  each  contain  two - f lash  1  amps  and  a 
7-mm-diameter-by-9-cm-long  Nd:YAG  rod  inside  a  BaS04-lined 
cavity.  The  oscillator  uses  a  telescopic  resonator  design  and  a 
flowing  liquid  dye  Q-switch,  as  shown  in  Figure  40.  The 
oscillator  output  is  coupled  to  the  amplifiers  by  the  polarizer. 
After  undergoing  a  single  pass  through  the  two  amplifier  heads, 
the  laser  pulse  is  reflected  by  a  phase  conjugate  mirror  (PCM) 
that  uses  SBS  in  high  pressure  CH4 .  The  SBS  "mirror"  generates  a 
retroref lected  wavefront-reversed  replica  of  the  incident  pulse 
with  a  reflectivity  of  ~80%  .  This  replica  passes  back  through 
the  amplifier  heads,  extracting  most  of  the  remaining  stored 
energy.  The  1.064-jUm  energy  output  is  ~700  mJ.  After  exiting 
through  the  polarizer,  the  output  is  frequency  doubled  using  type 
I  doubling  in  a  noncritical ly  phase  matched  CD*A  crystal.  Final 
output  is  ~250  mJ  at  532  nm . 

The  temporal  pulse  length  of  the  532  nm  pulse  is  ~19  ns  at 
FWHM ;  the  pulse  shape  is  shown  in  Figure  41.  The  beam  profile  is 
actually  somewhere  between  a  Gaussian  and  a  flat-top,  as  shown  in 
Figure  42,  a  plot  of  average  intensity  vs  radius  for  the  near- 
field  of  the  pump  beam.  The  pump  beam  was  found  to  be  SLM  in 
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Figure  39.  Low  energy  Nd : YAG  pump  laser  schematic. 

The  oscillator  operates  in  a  single¬ 
longitudinal  s ingle-transverse  mode. 
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Figure  40.  Single-longitudinal,  single-transverse 
mode  oscillator. 
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approximately  85%  of  its  pulses.  The  other  15%  of  the  pulses 
showed  modulation  of  the  temporal  shape,  indicating  the  presence 
of  two  or  more  modes.  The  pump  beam  linewidth  was  checked,  via  a 
measurement  of  its  coherence  length  (>3  m) ,  and  found  to  be 
<100  MHz.  This  implies  a  nearly  transform-limited  linewidth. 

The  experimental  arrangement  used  in  the  SLM  RA/SBS  experi¬ 
ments  is  shown  in  Figure  43.  Two  harmonic  beamsplitters  are  used 
as  the  first  two  turning  mirrors  to  eliminate  the  residual 
1.06  fim  energy  from  the  second  harmonic.  Next,  a  beamsplitter 
directs  ~30%  of  the  pump  energy  to  the  Raman  oscillator.  A  long¬ 
wave-pass  dichroic  beamsplitter  directs  this  pump  energy  through 
the  lens  into  the  Raman  oscillator.  The  dichroic  beamsplitter 
transmits  the  backward  1st  Stokes  Raman  energy  into  a  separate 
path  that  takes  it  through  the  Raman  amplifier.  After  its  first 
pass  through  the  amplifier,  where  it  picks  up  distortion  informa¬ 
tion,  the  backward  Stokes  "seed"  passes  through  another  dichroic 
beamsplitter  and  a  quarter-wave  plate  before  entering  the  PCM.  A 
phase  conjugate  replica  of  this  seed  proceeds  back  through  the 
amplifier  with  a  polarization  orthogonal  to  that  of  the  input. 
Because  of  the  conjugate  nature  of  the  returning  wave,  reciprocal 
aberrations  are  corrected  on  the  second  pass  through  the  ampli¬ 
fier.  On  its  return  path,  the  conjugate  beam  is  copropagating 
with  a  pump  beam  of  the  same  polarization. 

Thus,  the  Stokes  seed  experiences  gain  and,  possibly,  non- 
reciprocal  phase  aberrations.  After  the  second  pass,  the  ampli¬ 
fied  Stokes  is  extracted  from  the  original  beam  path  by  a 
polarizer  and  reaches  the  beam  diagnostics,  which  include  energy 
and  beam  profile  measurements.  Approximately  70%  of  the  original 
532  nm  pump  energy  traveled  through  a  delay  line  in  order  to 
reach  the  Raman  amplifier  at  the  same  time  the  Stokes  seed  was 
making  its  second  pass  through  the  amplifier.  Both  70%  of  the 
pump  energy  and  of  the  Raman  oscillator  outputs  are  imaged  to  the 
amplifier  and  reduced  in  size  in  order  to  increase  gain  and 
extraction.  This  imaging  of  their  respective  apertures  to  the 
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Figure  43. 


Experiment  used  to  demonstrate  nonreciprocal 
aberrations  using  the  single-mode  laser. 


,  f.  <r.  /,  <•. 


amplifier  eliminates  Fresnel  diffraction  until  the  beams  are 
propagating  through  the  amplifier.  Note  that  the  pump  and  Stokes 
seed  beam  polarizations  are  mutually  orthogonal  until  the  return 
pass  of  the  Stokes  seed  through  the  Raman  amplifier,  thereby 
minimizing  the  Raipan  gain  during  the  first  pass.  During  the 
experiments,  the  pump  and  seed  energies  were  varied  by  inserting 
neutral  density  filters  into  the  beams. 

The  backward  Stokes  Raman  oscillator  was  operated  with  an 
input  pump  energy  of  ~70  mJ  (~30X  threshold) .  When  seed  energy 
was  reduced,  the  neutral  density  filters  were  placed  in  the  red 
beam  only,  so  constant  pump  input  to  the  Raman  oscillator  was 
maintained.  Approximately  50%  energy  conversion,  or  65%  quantum 
conversion,  was  obtained.  The  temporal  pulse  shape  of  the  back- 


the  pump.  Nevertheless,  the  linewidth  is  less  than  one-third  of 
the  RA  gain  bandwidth.33  It  is,  however,  much  broader  than  the 
SBS  gain  bandwidth  in  methane34  (-20  MHz).  We  chose  acetone  as 
the  SBS  medium  for  these  experiments  because  it  has  an  SBS  gain 
bandwidth  of  -220  MHz.35 

4 . 5  EXPERIMENTAL  RESULTS 

Once  the  experimental  arrangement  was  assembled,  we  checked 
for  self-focusing  of  the  pump  beam  in  the  Raman  amplifier. 

During  these  experiments,  the  pump  beam  fidelity  transmitted  by 
the  RA  was  measured  as  a  function  of  pump  energy.  No  significant 
change  in  the  fidelity  was  observed.  The  pump  beam  quality  is 
given  in  Table  5. 

As  can  be  seen  from  the  table,  the  XDL  and  Strehl  calculated 
using  the  flat-top  beam  assumption  give  erroneous  results, 
because  our  pump  beam  has  a  smaller  diffraction  limit  than  a 


TABLE  5.  Characteristics  of  Single- 
Mode  Pump  Beam. 


PUMP  BEAM 
SINGLE-MODE 

GAUSSIAN 

BEAM 

FLAT-TOP 

BEAM 

NEAR  FIELD 
RADIUS  (cm) 

0.1430 

0.1572 

XDL 

1.58 

0.84 

STREHL 

0.54 

1.78 

perfect  flat-top  beam.  In  fact,  the  pump  beam  has  a  diffraction 
limit  only  1 . 5X  larger  than  a  Gaussian  beam.  In  what  follows,  w 
will  describe  RA/SBS  experiments  using  an  unaberrated  pump  beam, 
described  above,  and  an  aberrated  beam.  The  aberrator  used  for 
these  experiments  was  a  microscope  cover  slip  softened  by  flame 
and  allowed  to  distort  and  sag  over  a  piece  of  ceramic  material 
that  had  small  irregular  bumps  machined  into  it.  After  aberra¬ 
tion,  pump  beam  quality  decreased  to  a  Strehl  of  0.026,  or 
14.8  XDL. 
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STOKES 

UNABERRATED 

UNABERRATED 

ABERRATED 

ABERRATED 

SEED 

XDL 

STREHL 

XDL 

STREHL 

INPUT 

1.240 

0.781 

4.172 

0.150 

SEED  BEAM 

*0.034 

*0.029 

*0.034 

*0.002 

CONJUGATED 

1.681 

0.582 

2.471 

0.554 

SEED  BEAM 

*0.046 

±0.020 

±0.078 

*0.024 

From  the  table,  it  can  be  seen  that  the  phase  conjugate  fidelity 
was  not  perfect.  This  could  be  the  result  of  the  Stokes  seed 
linewidth  being  much  broader  than  the  SBS  gain  bandwidth.  The 
phase  conjugate  fidelity  of  the  narrower  bandwidth  pump  laser  was 
much  better  than  the  results  given  in  the  table  above. 

Having  established  the  quality  of  the  pump  beam  and  the  phase 
conjugation  fidelity,  we  proceeded  with  experiments  to  observe 
nonreciprocal  aberrations  in  the  RA/SBS  system.  In  Figure  45,  we 
plot  the  measured  gain  of  the  RA  for  both  aberrated  and  unaber¬ 
rated  beams.  For  an  aberrated  seed  beam  and  an  unaberrated  pump 
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Figure  45.  Raman  gain  as  a  function  of  the  pump  energy 

and  beam  aberration  using  the  single-mode  laser 
system.  Data  represented  by  the  triangles  and 
diamonds  are  for  an  aberrated  seed  beam;  the 
squares  and  the  diamonds  represent  the  aberrated 
pump  beam;  otherwise  the  beams  were  unaberrated. 


beam,  the  available  data  show  a  nearly  linear  increase  of  gain 
with  pump  energy.  This  indicates  that  saturation  of  the  pump  is 
playing  an  important  role  in  these  measurements.  A  least  squares 
fit  to  the  data  shows  that  (at  a  pump  energy  of  30  mJ)  gIL 
decreases  from  1.47  with  unaberrated  beams  to  0.91  with  aberrated 
beams.  In  Figure  46,  the  XDL  (Gaussian  beams  assumed)  of  the 
conjugated  seed  beam  is  plotted  for  aberrated  and  unaberrated 
beams.  As  expected,  the  XDL  of  the  amplified  seed  beam  increases 
with  pump  energy,  because  the  nonuniform  (Gaussian-1 ike)  pump 
beam  profile  gives  higher  gain  to  the  center  of  the  seed  beam, 
thereby  distorting  the  seed  beam  intensity  profile  and  increasing 
the  XDL.  When  the  pump  beam  is  aberrated,  the  unaberrated  seed 
beam  XDL  is  increased,  as  expected.  If,  however,  the  seed  beam 
is  aberrated,  we  find  a  reduction  in  the  XDL:  the  reason  is  not 
known . 

In  Figure  47,  we  plot  the  Strehl  corresponding  to  the  XDL  of 
Figure  46.  As  can  be  seen,  the  Strehl  decreases - monotonical ly 
with  increasing  pump  energy;  this  result  is  consistent  with 
Figure  46.  The  fact  that  the  Strehl  of  the  aberrated  and  unaber¬ 
rated  seed  beams  are  nearly  equal  is  a  result  of  our  aberration 
plate,  which  contains  primarily  short-scale-size  aberrations. 

The  resultant  aberrated  beam  consists  of  a  strong  central  peak 
with  a  broad  plateau  of  scattered  radiation.  Consequently,  the 
increase  in  XDL  is  much  larger  than  the  decrease  of  Strehl . 

In  the  single-mode  experiments  described  above,  the  RA  con¬ 
tained  a  pressure  of  468  psi  hydrogen,  while  the  Raman  oscillator 
(seed  generator)  contained  1252  psi.  This  differential  pressure 
was  chosen  to  exactly  compensate  the  SBS  shift  of  acetone 
(4.69  GHz). 36,37  For  comparison  with  the  simulations  described 
in  Section  2.3,  we  chose  to  vary  the  hydrogen  pressure  in  the 
Raman  oscillator.  We  chose  not  to  vary  the  RA  pressure  so  the 
Raman  gain  bandwidth  would  stay  fixed  at  ~1.5  GHz.  The  oscil¬ 
lator  pressure  was  varied  from  1099  to  1510  psi,  corresponding  to 
a  shift  of  -0.92  to  1.54  GHz  off  line-center.  The  measured  Raman 
gain  is  plotted  in  Figure  48.  The  data  show  an  unexpected 
resu 1 t . 
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Figure  46. 


Conjugate  beam  fidelity  as  a  function  of  pump 
energy  for  the  single-mode  RA/SBS  system.  The 
data  represented  by  the  triangles  and  diamonds 
are  for  an  aberrated  seed  beam;  the  squares  and 
the  diamonds  are  for  an  aberrated  pump  beam; 
otherwise  the  beams  were  unaberrated. 
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Conjugate  beam  fidelity  as  a  function  of 
pump  energy  for  the  single-mode  RA/SBS  system 
The  data  represented  by  the  triangles  and 
diamonds  are  for  an  aberrated  seed  beam;  the 
squares  and  the  diamonds  are  for  an  aberrated 
pump  beam;  otherwise  the  beams  were 
unaberrated . 
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Figure  48.  Raman  gain  as  a  function  of  the  RA  cell 

pressure  using  the  single-mode  RA/SBS  system 
The  beams  are  unaberrated  except  for  the 
squares,  where  the  seed  beam  was  aberrated. 
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The  measured  Raman  gain  varied  linearly  with  the  oscillator 
pressure.  The  expected  result  was  that  the  gain  would  peak  on 
"line  center,"  or  at  about  1252  psi .  Our  results  indicate  that 
the  maximum  Raman  gain  occurs  with  an  oscillator  pressure  of 
~1500  psi.  This  corresponds  to  a  SBS  shift  ~1.3X  larger  than 
expected,  or  a  pressure  dependence  of  the  Raman  shift  ~0.76X 
smaller  than  expected.  The  interplay  of  the  pump  beam,  seed 
beam,  and  RA  linewidths  could  also  contribute  to  these  unexpected 
results.  Because  of  time  constraints,  we  were  unable  to  pursue 
this  issue  further.  We  did,  however,  measure  the  conjugate  beam 
fidelity  as  a  function  of  the  Raman  oscillator  pressure.  These 
results  are  shown  in  Figures  49  and  50.  In  Figure  49,  the  con¬ 
jugate  beam  XDL  is  plotted  as  a  function  of  oscillator  pressure, 
and  in  Figure  50  the  conjugate  Strehl  is  plotted.  No  significant 
change  in  conjugate  beam  fidelity  is  observed  as  the  Raman  oscil¬ 
lator  pressure  is  changed. 

In  summary,  we  have  assembled  a  RA/SBS  system  to  evaluate 
nonreciprocal  aberrations  using  a  single-longitudinal/single- 
transverse-mode  laser.  The  diffraction  limit  of  the  RA/SBS 
system  was  significantly  better  than  is  theoretically  possible 
for  flat-top  beams,  and  was  only  ~1.5X  diffraction  limited  using 
Gaussian  beam  calculations.  Using  fixed  aberrators,  we  demon¬ 
strated  that  SBS  phase  conjugation  was  capable  of  good  compensa¬ 
tion  for  large  induced  reciprocal  aberrations.  The  conjugate 
beam  fidelity  was  then  measured  as  a  function  of  the  RA  gain  for 
both  aberrated  and  unaberrated  beams.  We  found  that  the  ampli¬ 
fied  seed  beam  fidelity  decreased  monotonically  with  increasing 
pump  energy,  probably  as  a  result  of  nonuniform  gain  resulting 
from  the  Gauss i an-1 ike  intensity  profile  of  the  pump  beam. 
Investigation  of  the  conjugate  beam  fidelity  as  a  function  of 
frequency  mismatch  between  the  conjugate  seed  and  the  peak  of  the 
Raman  gain  yielded  unexpected  results:  further  investigation  is 
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Conjugate  beam  fidelity  as  a  function  of 
Raman  cell  pressure  for  the  single-mode  RA/SBS 
system.  The  circles  represent  unaberrated  seed 
beam  results  and  the  squares  aberrated  seed 
beam  results.  The  pump  beam  was  unaberrated  in 
both  cases . 


SECTION  5 


SUMMARY  AND  RECOMMENDATIONS  FOR  FUTURE  WORK 

5.1  COMPARISON  OF  SIMULATIONS  WITH  SINGLE-MODE  LASER  RESULTS 

As  discussed  in  Section  2,  we  have  simulated  the  conjugate 
beam  fidelity  for  a  situation  that  closely  corresponds  to  the 
single-mode  laser  data  discussed  in  Section  4.4.  These  results 
are  shown  in  Figures  51  and  52.  In  Figure  51,  the  calculated 
reduction  of  Strehl  is  plotted  as  a  function  of  the  small-signal 
gain  parameter.  The  solid  lines  show  the  simulation  results, 
while  the  open  circles  and  squares  show  single-mode  experimental 
results  for  an  aberrated  and  unaberrated  seed  beam,  respectively. 
The  pump  beam  was  unaberrated  in  both  cases .  As  were  the 
simulation  results,  the  experimental  results  are  normalized  to 
the  Strehl  of  the  unamplified  seed  beam.  The  figure  shows  that 
our  experimental  results  lie  between  the  simulation  of  a  perfect 
conjugator  and  the  phase-only  conjugator.  This  result  implies 
that  the  fidelity  of  our  SBS  phase  conjugator  is  not  perfect,  but 
also  that  the  impact  of  a  nonreciprocal  aberration  (nonuniform 
gain)  is  as  strong  as  the  simulations  predict.  In  Figure  52,  the 
experimental  results  are  plotted  along  with  the  simulations  of 
RIVED.  As  seen  in  the  figure,  all  our  experimental  results  lie 
at  pump  energies  too  small  for  the  RIVED  effect  to  be 
significant.  In  fact,  the  simulation  results  correspond  to  the 
reduction  of  Strehl  at  the  end  of  the  pulse,  the  effect  being 
absent  at  the  beginning.  As  such,  the  effective  energy  for  the 
experimental  results  should  be  reduced  by  about  a  factor  of  2 
(assuming  linear  reduction  of  Strehl  with  cumulative  energy 
during  the  pulse) . 

Comparison  of  the  multimode  laser  results  with  the 
simulations  was  not  done  because  the  simulation  does  not  take 
into  account  the  broad  laser  bandwidth. 
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Figure  51.  Comparison  of  simulated  and  experimental  conjugate 
beam  fidelity  as  a  function  of  RA  gain.  Simulation 
results  are  represented  by  the  solid  symbols:  the 
solid  dots  are  for  an  unaberrated  seed  beam.  The 
solid  diamonds,  squares,  and  triangles  are  for  an 
aberrated  seed  beam  with  perfect  conjugation,  phase 
only  conjugation,  and  without  conjugation, 
respectively.  In  all  cases,  the  pump  beam  was 
assumed  to  be  an  unaberrated  Gaussian.  The  open 
symbols  show  the  single-mode  experimental  results, 
normalized  to  the  seed  fidelity,  as  were  the 
simulations.  The  open  circles  and  squares  are  for 
the  aberrated  and  unaberrated  seed  beam, 
respectively . 
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5.2  RECOMMENDATIONS  FOR  FUTURE  WORK 

During  this  program,  we  demonstrated  the  reduction  of  beam 
fidelity  with  Raman  gain  in  a  RA/SBS  system.  These  reductions 
result  from  nonreciprocal  aberrations  in  the  RA.  Theoretical, 
numerical  simulations  and  experimental  results  are  in  good 
qualitative  agreement.  Quantitative  agreement  between  single¬ 
mode  experiments  and  the  simulations  is  within  experimental 
uncertainties.  The  RA/SBS  system  shows  great  promise  for  high 
energy  laser  applications,  but  its  implementation  requires  a 
deeper  understanding  of  the  basic  device  physics.  Within  the 
scope  of  this  program  we  were  able  to  demonstrate  the  concept, 
and  the  impact  of  nonreciprocal  aberrations,  and  to  identify 
nonreciprocal  processes.  The  experimental  implementation 
revealed  knowledge  gaps  for  the  basic  processes  in  the  nonlinear 
devices.  These  knowledge  gaps  reveal  clear  avenues  for  future 
investigation : 

•  SBS  phase  conjugators: 

Conjugation  fidelity  for  wide  bandwidth  pump  beams 
Conjugation  fidelity  for  highly  aberrated  pump  beams 
Phase  vs  amplitude  conjugation 

•  Backward  stimulated  Raman  oscillators: 

Linewidth  of  generated  seed 
Dependence  on  pump  laser  bandwidth 
Dependence  on  aberration  on  pump  beam 

•  Raman  Amplifier: 

Reciprocal  aberrations  induced  by  thermal  effects 
Operation  in  high  Fresnel  number  geometry. 

In  addition,  simulations  of  transient  effects,  such  as  RIVED,  are 
needed  to  correctly  determine  the  effect  of  this  aberration  in 
pulsed  laser  systems. 
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APPENDIX  A 


PROPAGATION  MODELING 


The  Basic  Approach 

The  propagation  phenomena  associated  with  Raman  scattering 
can  be  modeled  with  essentially  the  same  techniques  we  have  used 
to  model  the  effects  of  turbulence  and  thermal  blooming  for 
propagation  in  the  atmosphere  and  fields  within  laser  resonators. 
The  basic  idea  is  that  the  effects  of  the  medium  can  be 
approximated  by  a  series  of  localized  layers  of  gain  and  phase 
shift.  The  wave  is  propagated  from  one  layer  to  the  next  with  a 
free-space  propagator  that  solves  the  wave  equation  in  the 
paraxial  approximation,  which  is  done  with  a  Fast  Fourier 
Transform  (FFT)  algorithm.  This  approach  is  highly  accurate  if 
the  distance  between  the  layers  is  not  too  large  and  the  wave  is 
adequately  sampled  in  the  direction  transverse  to  the 
propagation. 

The  principal  difference  between  the  requirements  for 
modeling  Raman  scattering  phenomena  and  for  propagation  modeling 
is  the  necessity  in  the  Raman  problem  to  include  the  effect  of 
more  than  one  field.  Where  all  the  fields  are  propagating  in 
nominally  the  same  direction  (e.g.,  in  an  amplifier  in  which  the 
effects  of  backward  waves  are  neglected),  the  necessity  to 
include  more  than  one  field  impacts  the  modeling  only  by  an 
approximately  linear  increase  in  computational  time  and  memory 
requirements.  No  significant  changes  in  the  modeling  technique 
are  required.  Where  counterpropagating  fields  exist  (e.g.,  in  an 
amplifier  where  backward  wave  scattering  exists),  the  difficulty 
of  the  modeling  is  significantly  increased  by  virtue  of  the  fact 
that  the  boundary  conditions  on  the  forward  and  backward  waves 
are  not  given  at  the  same  location.  The  forward  waves  are 
specified  at  one  end  of  the  active  medium  and  the  backward  waves 
at  the  other.  We  have  previously  explored  various  techniques  for 
treating  such  "two-point"  boundary  value  problems  and  are 
actively  investigating  new  techniques  under  IR&D  support. 


The  propagation  from  one  gain/phase  layer  to  the  next  is 
calculated  using  a  converging  coordinate  system.  This  allows  us 
to  consider  cases  in  which  the  fields  are  strongly  converging  or 
diverging  as  well  as  collimated.  The  converging  coordinate 
system  is  defined  in  terms  of  the  convergence  (or  divergence)  of 
a  Gaussian  beam  having  an  e-folding  radius  p0  equal  to  one-third 
the  diameter  of  the  fields  being  propagated.  With  this  choice  of 
e-folding  radius,  the  fictitious  Gaussian  beam  has  approximately 
90%  of  its  energy  confined  within  the  diameter  of  the  propagated 
field.  The  converging  coordinate  variables  are 


P 0p(z) 


9  oP(z) 


P(z)  =  a'1'2  [(1  -  z/f)2  +  a2(z/kp2)2] 1/2 


(A.  1) 


(A.  2) 


(A- 3) 


pTTT 


(A.  4) 


=  *  (f /kp*) 


(A. 5) 


(A.  6) 


where  x,y,z  are  Cartesian  coordinates  with  z  the  direction  of 
propagation,  f  the  distance  to  the  focus  of  the  converging 
Gaussian  beam,  a  a  parameter  that  allows  us  to  consider  cases  in 
which  the  size  of  the  focal  spot  may  exceed  the  diffraction  limit 
(e.g.,  where  the  initial  fields  are  distorted),  and  d  the 
diameter  of  the  beams  in  the  plane  z  -  0. 


When  expressed  in  terms  of  the  converging  coordinate  system 
defined  above,  the  free-space  propagator  satisfies  the  equation 


du 


0 


(A  .  7) 


In  most  applications,  the  increment  in  £  in  going  from  one 
gain/phase  layer  to  the  next  is  small.  Thus,  we  approximate  the 
solution  of  Eq.  (A. 7)  by 

u  2  exp  [- i  +  d)  €/2J  ^  (A. 8) 


where  v  satisfies 

di/  i  (d*v  a2M  _  n 

d(  2  +  a^J  '  0 


(A.  9) 


Equation  (A. 9)  can  be  easily  solved  with  a  FFT  routine.  The 
solution  in  the  transform  domain  is 


V(K,£  +  A£)  =  V(K.O  exp  ( - iK2 A£/2)  ,  (A. 10) 

where  V(K,£)  is  the  Fourier  transform  of  the  field  at  the  plane 
£.  In  cases  where  the  distance  between  the  gain/phase  layers  is 
larger,  we  solve  Eq .  (A. 7)  directly  by  expanding  ui  in  Hermite 

polynomials.  We  refer  to  this  solution  as  a  "Hermite 
propagator,"  and  to  the  FFT  solution  as  a  "FFT  propagator."  In 
most  cases  we  use  the  FFT  propagator  because  it  is  more  efficient 
computationally.  The  Hermite  propagator  is  sometimes  used  in 
optical  train  calculations  and  in  atmospheric  propagation 
calculations  in  which  the  wave  is  propagated  in  one  large  step 
from  the  top  of  the  atmospheric  layer  to  a  distant  point  in 
space.  Although  the  Hermite  propagator  is  available  if  needed, 
we  anticipate  that  the  FFT  propagator  will  be  used  almost 
exclusively  in  the  Raman  modeling  studies. 


Sampling  and  Storage  Requirements 


The  number  of  mesh  points  required  in  the  propagation  mesh 
depends  on  the  highest  significant  spatial  frequency  of  the 
variation  of  the  complex  field  amplitude  in  the  direction 
transverse  to  the  propagation.  Because  we  were  interested  in  the 
effects  associated  with  dirty  pump  waves,  it  was  necessary  to  use 
a  mesh  having  at  least  128  sample  points  for  each  transverse 
direction  (i.e.,  a  128  by  128  mesh  for  a  3-D  (x,y,z)  model  and  a 
128  points  mesh  for  a  2-D  (x,z)  model) .  Previous  experience  with 
the  modeling  of  the  distorted  output  of  real  high  power  laser 
devices  has  indicated  that  more  sample  points  could  be  needed  if 
the  code  is  to  model  the  effects  of  heavily  distorted  pump  waves. 

Having  selected  the  number  of  points  in  the  propagation  mesh, 
we  must  consider  the  impact  on  the  requirements  for  the  number  of 
steps  to  be  taken  along  the  direction  of  propagation.  In  order 
to  adequately  sample  the  fields  in  the  transform  domain,  we  must 
take  enough  propagation  steps  that  the  phase  factor  in  Eq .  (A. 10) 

changes  by  no  more  than  a  radian  or  so  between  adjacent  mesh 
points.  The  spatial  wavenumber  K  is  related  to  the  number  of 
mesh  points  N  and  the  increment  A£  in  the  transverse  variable  £ 


(K„,  K  )  =  [ (m-N/2) AK ,  (n-N/2)AK] 


AK  =  27r/NA$ 


(A. 11) 


(A. 12) 


It  can  readily  be  shown  that  the  maximum  phase  difference  between 
adjacent  mesh  points  is 


27T2  A< 


N  (A$) 


(A. 13) 


It  is  immediately  obvious  from  this  relation  that  one  must 
decrease  the  distance  between  the  gain/phase  layers  when  the 
distance  between  the  transverse  mesh  points  is  decreased.  The 


WTO 


increments  A£,  A£  are  related  to  the  number  of  propagation  steps 
N*t*p  anc^  tbe  number  of  mesh  points  Nmp  distributed  across  the 
diameter  of  the  input  beams  by  the  relations  (note  that  Nmp<N) 


AS 


3a1' 2 
Nmp~  1 


(A. 14) 


N.t.p 


(A. 15) 


For  a  case  in  which  the  active  medium  is  2  m  long  and  the 
input  waves  are  collimated  (f  =  «) ,  with  d  =  0.5  cm,  X  =  0.4  pm, 
N  =  128,  Nmp  =  89,  and  a  -  1,  the  maximum  phase  shift  between 
mesh  points  is 


A4>m,x  =  0.05  N2p/NNstep  =  3/Nstep  .  (A. 16) 

Thus,  something  of  the  order  of  ten  propagation  steps  would  be 
required  to  yield  adequate  sampling  in  the  transform  domain  in 
this  case.  Of  course,  we  must  also  consider  the  number  of 
gain/phase  layers  required  to  approximate  the  continuous 
distribution  of  gain  (phase  shift)  introduced  in  the  real  medium. 
Generally,  the  number  of  gain/phase  layers  (which  is  equal  to 
N.t.p+1)  must  be  sufficient  to  limit  the  gain  introduced  by  a 
single  layer  to  unity  or  less  . 


The  storage  requirements  are,  of  course,  intimately  tied  to  the 
sampling  requirements.  They  also  depend  on  whether  the  fields 
are  all  propagating  in  the  same  direction  or  not.  The  amount  of 
memory  available  on  a  Cray  is  adequate  to  handle  most  cases  of 
interest . 


Distorted  Wave  Modeling 


We  have  previously  developed  techniques  for  generating  2-D  random 
distributions  having  specified  spatial  correlation  properties. 
These  techniques  are  used  extensively  in  our  adaptive  optics  and 


s 


propagation  software  to  model  the  effects  of  atmospheric 
turbulence .  The  same  techniques  can  be  used  to  generate 
realizations  of  distorted  waves  in  the  Raman  modeling  studies. 

For  example,  "dirty"  pump  waves  and  replicas  of  atmospherically 
degraded  waves  for  use  in  the  Raman  amplifier  studies  can  be 
generated  with  these  techniques.  Generally,  a  realization  of  a 
2-D  random  distribution  H(x)  having  the  spatial  spectrum  <p(K)  can 
be  obtained  from  the  relation 

H(x)  =  /  dS(K)  exp(iKx)  (A. 17) 


with 


<dS (K) >  =  0. 


(A. 18) 


V. 

>. 


<dS(K1)dS(K2)>  =  0. 


Ki  t  K2 


(A. 19) 


=  dKj4>  (Kj )  Kx  =  K2 


V  4 


In  atmospheric  propagation  studies  we  use  the  following  modified 
Kolmogorov  spectrum  to  describe  the  phase  fluctuations  associated 
with  turbulence: 


<J>  (K)  =  0.207C2k2AzK0-  1  1  /3  (1  +  K2/K2)'11/6 


(A. 20) 


We  have  also  used  the  same  technique  to  generate  realizations  of 
waves  having  amplitude  scintillations  associated  with  atmospheric 
turbulence.  Dirty  pump  waves  can  be  modeled  by  similar 
techniques  in  which  the  amplitude  and  phase  variations  of  these 
waves  are  described  by  appropriate  spatial  spectra. 
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APPENDIX  B 


PHASE  ERRORS  IN  RAMAN  AMPLIFIERS 

As  an  inherent  part  of  the  Raman  amplification  process, 
energy  is  deposited  in  the  Raman  medium  in  the  form  of 
vibrational  excitation.  This  vibrational  energy  thermalizes  and 
appears  as  heat.  As  is  well  known,  a  homogeneously  heated  gas 
can  induce  phase  aberrations  in  an  optical  beam.  In  Raman 
amplifier  design,  a  technique  must  be  found  to  eliminate  these 
thermally  induced  aberrations  in  order  to  permit  PRF  scaling 
while  maintaining  acceptable  beam  quality  and  pointing  stability. 
In  this  appendix,  we  estimate  the  magnitude  of  the  thermal  phase 
distortion . 

Consider,  as  a  typical  case,  a  hydrogen  Raman  amplifier 
pumped  by  a  5-J,  10-Hz,  525-nm  beam.  The  first  vibrational  state 
of  hydrogen  has  an  energy  ~20%  of  a  green  photon;  hence  the 
energy  deposited  in  the  gas  is  approximately  20%  of  the  pulse 
energy,  or  ~1  J.  Operation  at  10  Hz  corresponds  to  an  average 
heat  load  in  the  gas  of  ~10  W.  This  heating  induces  temperature 
gradients  that  lead  to  density  fluctuations  through  a  variety  of 
mechanisms,  such  as  thermal  expansion  and  free  convection. 

Temperature  variations  in  the  Raman  medium  may  be  turbulent 
or  laminar.  Laminar  variations  are  typically  larger  in 
magnitude,  and  of  larger  spatial  scale,  while  the  turbulent 
fluctuations  are  smaller,  with  a  small  spatial  scale.  The  phase 
distortions  from  laminar  temperature  variations  are  typically 
tilt,  focus  (or  defocus),  and  astigmatism,  while  the  turbulent 
distortions  are  small  scale  random  fluctuations  that  can  cause 
wide  angle  scattering  when  the  beam  is  propagated  to  the  far 
field. 

Detailed  thermal  calculations  for  a  typical  situation  have 
been  performed.  Assuming  a  2-cm-d i ameter ,  50-cm-long  column  of 
the  Raman  gas  at  60  atm  pressure  is  uniformly  heated  with  10  W, 
we  find  that  the  temperature  profile  changes  from  RMS 


fluctuations  of  a  few  tenths  of  a  degree  in  a  medium  with  forced 
turbulence  (i.e.,  with  an  internal  fan),  to  a  roughly  parabolic 
profile  with  a  peak  temperature  change  of  about  10°  for  a  laminar 
medium . 

Temperature  changes  alter  the  local  gas  density.  According 
to  the  ideal  gas  law,  AN/N  =  -AT/T  where  N  is  the  molecular 
number  density  and  T  is  the  temperature.  The  refractive  index 
change  accompanying  this  density  change  is  found  using  the 
Lorentz-Lorenz  equation 


n2  -  1 
n2  +3 


4ttN 


(B.l) 


where  n  is  the  refractive  index  and  Q  is  the  molecular 
polarizability.  This  equation  can  be  expressed  in  a  more 
convenient  form 


n 


-  1  = 


4?rNa 


1  -  (4tt/3)N  a 


(B  •  2) 


where  a  pressure  of  60  atm  is  assumed. 


The  RMS  phase  error  a  resulting  from  the  An  is  given  by 

a2  ~  ck2  L  L0  <An2>  ,  (B.7) 

where  c  is  a  constant  ~1,  k  =  2ir/X  ,  X  is  the  optical  wavelength, 
L  is  the  length  of  the  medium,  L0  is  the  outer  scale  of  the 
turbulence,  and  <An2>  is  the  variance  of  the  refractive  index. 

This  expression  can  be  derived  by  considering  the  random 
phase  aberration  A$  due  to  a  random  refractive  index  fluctuation 
An 


A0  =  k  Jgdz  An  ,  (B.8) 

and  finding  the  RMS  a  by 

o 2  =  <A <p2>  =  k2  c/^dzAnCz)  /kdz’An(z’)>  ,  (B.9) 

where  the  brackets  denote  a  time  average.  If  the  refractive 
index  fluctuations  have  an  exponential  correlation  function  of 
the  form 

<An(z)  An(z’)>  =  <An2>  exp (- j z-z 0 | /L0 )  ,  (B.10) 


then 


a2  =  2k2LL0  <An2>  (B . 1 1 ) 

Other  correlation  functions  give  the  same  functional 
dependence  and  numerical  coefficients  ranging  from  x/4  to  ir/2.38 
The  phase  error  can  now  be  calculated  by  substituting  Eqs .  (B.5) 

38L.A.  Chernov,  Wave  Propagation  in  a  Random  Medium  (Dover,  NY, 
1967) ,  75. 


and  (B.6)  into  (B.ll),  and  using  L  =  50  cm,  L0  =  0.3  cm  (the 
outer  scale  is  typically  ~0 . IX  the  tube  diameter,  assumed  to  be 
~2  cm).  Considering  first  the  laminar  case,  we  calculated 
AT~10*K ,  and 

a  =  150  rad  (hydrogen)  (B.12 

a  =  470  rad  (methane)  (B.13 

For  the  turbulent  case,  AT  was  found  to  be  ~0.1®K,  and  the 
resulting  phase  errors  are  smaller  than  the  above  values  by 
factor  of  100. 
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